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I Foreword 

This  report  Is  based  on  a lecture  series  presented  by  the 
authors  at  the  von  Karma'll  Institute  for  Fluid  Dynamics,  Rhode- 
Saint-Genese,  Belgium, >in  May  1965;  and  at  the  University  of 
Maryland,  College  Park,  Maryland,  ln^July  1965.  The  lectures  were 
prepared  under  the  Joint  auspices  of  the  David  Taylor  Model  Basin 
and  the  Naval  Air  Development  Center.  They  were  presented  In 
Belgium  under  the  joint  sponsorship  of  the  von  Karman  Institute 
and  the  Advisory  Group  for  Aerospace  Research  and  Development 
(AGARD) ; and  in  Maryland  under  the  sponsorship  of  the  Assistant 
Secretary  of  the  Navy  for  Research  and  Development. 

) * The  revised  lectures  will  be  presented  as  follows: 

> • (±) 

A.  Introductory  Survey 

B.  'Air  Cushion  Mechanics  > (3  ) ^ 

C.  Internal  Aerodynamics.  ^ ?,) 

D.  Drag  / 


F.  ^Cushion  Contributions  to  Stability 

G.  Seakeeping  , 

H.  Performance  Summary 
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NOTATION 

cushion  length,  ft 

cushion  beam,  ft 

cushion  area,  ft3 

cushion  perimeter,  ft  * 

"daylight"  gap,  ft 

cushion  pressure,  lb/ ft3  (gage) 

total  pressure  at  nozzle  exit,  lb/ft3  (gage) 

gross  weight,  lb 
gap  area  (hC) , ft3 

forward  velocity,  ft/sec 


forward  velocity,  knots 
cushion-reference  velocity 

air  density,  slugs/ft3 

flow  quantity  rate,  ft3/sec 
total  power,  lb-ft/sec 
propulsion  power,  lb-ft/sec 

cushion  power,  lb-ft/sec 


ideal  total  power,  lb-ft/sec  (corresponding  to  efficien- 
cies of  unity) 

ideal  propulsion  power,  lb-ft/sec 


ideal  cushion  power,  lb-ft/sec 


dynamic  pressure  (pVQ‘?/2),  lb/ ft* 

/ D , 


W 


total  drag  coefficient 
total  drag,  lb 

coefficient  of  "other"  drag  excluding  wave  drag 
wavemaking  drag,  lb 


fD  - Dw 
V qaS 


discharge  coefficient 


(Q  + s8  vc) 


ill 
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INTRODUCTION 


It  Is  conventional  to  begin  any  discussion  of  Ground  Effect 
Machines  with  a reminder  that  these  vehicles  have,  as  yet,  enjoyed 
only  a very  short  history  of  development.  While  this  might  become 
tiresome  to  those  who  interest  themselves  in  very  many  3uch  discus- 
sions, it  is  nevertheless  one  of  the  facts  which  it  is  absolutely 
essential  to  keep1  in  mind,  if  one  is  to  acquire  or  maintain  a rational 
perspective  of  this  important  new  field. 

The  history  of  serious  study  and  development,  having  pertinence 
to  the  current  state  of  affairs,  goes  back,  in  fact,  only  about  eight 
or  nine  years;  to  Cockerell's  initial  success  in  attracting  official 
support  for  his  ideas,  in  the  United  Kingdom,  and  to  the  beginnings  of 
serious  research  programs  in  the  United  States.  From  a purely  histor- 
ical point  of  view,  of  course,  one  must  go  back  much  further,  perhaps 
a century  or  more,  to  find  the  first  of  a remarkable  number  of  inde- 
pendent conceptions  of  GEM's  of  various  types;  but  these  seem  to  have 
had  little  or  no  influence  on  the  modern  development  of  the  field. 

It  is  not  very  meaningful,  of  course,  to  compare  the  current 
state  of  GEM  development  with  that  of  the  airplane  in  the  year  1911, 
as  is  sometimes  done.  The  environment  for  technological  development, 
and  the  pace  at  which  it  is  possible  to  carry  such  developments  forward, 
have  changed  a very  great  deal  since  the  Wright  brothers  flew  their 
first  airplane.  Even  so,  nine  years  (or  more  appropriately,  six  years 
since  SR-Nl  made  her  first  flight)  is  still  a very  short  time,  indeed, 
in  the  life  of  an  important  new  form  of  transportation.  Some  conse- 
quences of  this  youthfulness  of  the  field  are: 

1.  The  state  of  the  art  must  be  regarded  as  primitive  (the  remark- 
able practical  success  of  some  of  the  British  Hovercraft  notwithstanding). 

2.  Conclusions  and  opinions  of  almost  every  kind  should  be  reached 
and  held  tentatively.  Physical  principles  will  not  change,  of  course; 
but  the  apparent  practical  significance  of  some  of  che  most  fundamental 
principles  has  been  changing  rapidly,  and  this  must  be  expected  to  con- 
tinue as  the  field  continues  to  develop  and  expand. 


I 


3.  In  common  with  moat  fields  of  technology,  whether  primitive 
or  advanced,  the  most  important  facts  known  about  Ground  Effect  Ma- 
chines are  rather  simple,  elementary  facts;  but  also,  the  most  important 

of  the  questions  which  remain  unanswered  are  rather  simple,  elementary 

I 

[questions.  It  follows  that  the  role  of  subtle  and  sophisticated  mathe- 
matical analysis  is,  as  yet,  comparatively  unimportant  in  this  field. 

\ The  simple,  commonsense  approach  is  usually  far  more  rewarding;  and 

even  the  simplest  approaches  lead  quickly  into  areas  where  the  uncer- 
tainties of  the  underlying  experimental  facts  exceed  the  inherent  mathe- 
matical errors  of  first-order  analysis. 

We  will  try  to  follow  the  simplest  and  clearest  paths  available, 
seeking  to  expose  the  main  ideas  in  their  elementary  forms,  and  relying 
on  intuitive  plausibility  in  lieu  of  proof. 

In  fact,  if  we  wish  to  pursue  our  subject  to  the  point  of  conclu- 
sions, we  have  little  choice  but  to  rely  on  intuitive  approaches  at 
many  stages  along  the  way.  But  let  the  reader  beware!  Conclusions 
reached  in  such  a manner  are  not  to  be  accepted  as  permanent  truths. 

They  are  to  be  regarded  as  a form  of  tentative  interpretation  of  available 
facts,  which  must  be  constantly  reinterpreted  as  new  facts  come  along. 

SCOPE 

Not  the  least  of  the  vagaries  of  the  Ground  Effect  Machine  field 
i3  the  meaning  of  the  term  Ground  Effect  Machine  (GEM).  This  term  is, 
in  fact,  only  one  of  a number  of  terms  which  have  been  used  rather 
loosely  and  interchangeably  in  reference  to  the  same  general  class  of 

vehicles.  A partial  listing  of  some  of  the  more  common  terms  is: 

* 

Ground  Effect  Machine 

Hovercraft 

Air  Cushion  Vehicle 

Surface  Effect  Ship 

Hydroskimmer 

Various  authors  have  elected  to  consider  one  or  another  of  these  terms 
to  be  descriptive  of  a very  broad  class  of  vehicles,  using  others  to 
describe  subclasses  within  the  broad  class;  but  there  seems  to  have  been 
little  or  no  agreement  on  this  among  different  authors. 


; 


We  elect  to  consider  that  all  of  the  above  terras  are  synonymous, 
describing  the  same  broad  class  of  vehicles,  which  we  proceed  to  define 
as  follows  : 

Ground  Effect  Machine  (GEM)  — Any  vehicle  which  causes  air  pressure 
higher  than  atmospheric  to  be  exerted  on  the  surface  of  the  earth  di- 
rectly below  it,  so  that  the  integral  of  this  additional  air  pressure, 
over  the  area  of  the  vehicle  planform  projected  on  the  earth's  surface, 
produces  a force  comparable  to  (say,  greater  than  half)  the  total  weight 
of  the  vehicle.  The  air-filled  space  between  the  earth's  surface  and 
the  vehicle  is  called  the  "air  cushion." 

The  definition  embraces  a much  broader  class  of  vehicles  than  we 
will  consider  in  detail.  We  will,  in  fact,  consider  in  detail  only 
vehicles  with  air  cushions  of  the  plenum  and  peripheral  jet  types,  and 
certain  design  variations  of  these  involving  sidewalls  and  other  forms 
of  mechanical  "cushion  containers.”  (This  narrower  class  of  vehicles 
includes,  for  example,  the  Westland  Hovercraft,  3R.N1  through  SR.N6;  the 
Denny  sidewall  craft  D-l  and  D-2;  the  U.  S.  Navy's  XR-1  "Captured  Air 
Bubble";  the  Bertin  B—  6 ; and  a large  number  of  similar  craft.) 

Before  embarking  on  detailed  consideration  of  the  fluid  dynamic 
design  principles,  however,  it  will  be  profitable  to  consider  briefly 
some  of  the  more  obvious  characteristics,  similarities,  and  differences 
of  various  types  of  craft,  including  some  which  lie  beyond  the  scope  of 
what  we  will  later  consider  in  detail. 

The  plenum  type  craft  is  comparatively  very  easy  to  understand, 
yet  possesses  most  of  those  characteristics  of  Ground  Effect  Machines 
which  must  be  understood,  if  one  is  to  understand  Ground  Effect  Ma- 
chines at  all.  We  will  therefore  consider  the  Plenum  first,  and  in 
somewhat  more  detail  than  other  types  of  craft,  for  purposes  of  this 
introductory  survey. 

SIMPLE  PIJifTUM  CRAFT 


Consider  a craft  as  sketched  below: 


The  cushion  is  contained  in  a large  cavity,  between  the  ground,  the 
base  of  the  craft,  and  the  "plenum  walls"  extending  down  from  the  base. 
Air  is  pumped  directly  into  the  cushion  through  a compressor,  and 
escapes  through  a "daylight"  gap  at  the  lower  edge  of  the  cushion,  at 


velocity  V£  . If  we  choose  the  primary  variables  to 

be: 

Cushion  Area 

s 

ft2 

Cushion  Perimeter 

C 

ft 

Daylight  Gap 

h 

ft 

Cushion  Pressure 

Ap 

psfg 

Discharge  Coefficient 

£> 

c 

— 

Other  properties  of  the  system  follow 

directly 

Lift  = Weight 

Wifip  S 

lb 

Gap  Area 

S - hC 
g 

ft2 

Reference  Velocity 

V = /-  Ap 

c J Pa 

ft/sec 

Flow  Quantity 

Q = V S £ 
c g c 

ft3 /sec 

Ideal  Cushion  Power  ] 

P = Ap  Q 

lb-ft/sec 

J 
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These  relationships  become  slightly  modified  with  forward  speed, 
Vq  , when  the  dynamic  pressure 


p V s/2 
Ka  o 


becomes  comparable  in  magnitude  to  the  cushion  pressure  Ap.  At  mod- 
erate speeds,  however,  we  can  write: 

Cushion  Specific  Power  (Ideal) 


P . S 

ci±  _ _£  * 2. 

W V Sc' 

o 


J qfl/Ap  ^Plenum,  _a  < < 1^  [A 


The  ideal  propulsion  power  is  the  product  of  total  drag  and 
velocity.  If  we  define  a total  drag  coefficient 


D q S 
na 


then 


P . = q S V 

p,i  Da  o 


and  we  can  write 

Propulsion  Specific  Power  (Ideal) 


JLil  = £ = c — 
W V W D Ap 

o 


[A 


also, 


Total  Specific  Power  (Ideal) 


I 
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The  specific  power  (often  called  "equivalent  drag/lift  ratio")  la 
one  of  the  most  important  of  the  various  nondimens lonal  vehicle  perform- 
ance parameters.  For  a vehicle  of  given  gross  weight  W,  the  weight  of 

fuel  W,  consumed  during  a trip  of  a given  distance  is  directly  propor- 

f P 

tional  to  the  total  specific  power  (see  Figure  A-l). 

o 

The  simple  form  of  Equations  [A-l]  through  [A-3]  conceals  many 
complications  with  which  we  will  deal  in  the  succeeding  sections.  The 
cushion  pressure  and  cushion  power  can  vary  somewhat  with  forward 
speed;  and  the  total  drag  coefficient  is  composed  of  many  components 
(external  aerodynamic  drag,  momentum  drag,  wavemaking  drag,  hydrodynamic 
drag,  etc.)  which  vary  in  different  ways  with  speed.  (Note,  is 
constant  if  D is  proportional  to  VQ2.)  Also,  we  have  only  the  ideal 
specific  power,  with  no  allowance  yet  for  the  efficiency  of  the  pro- 
peller, compressor,  ducts,  etc. 

Nevertheless,  an  excellent  beginning  toward  insight  into  the  nature 
of  Ground  Effect  Machines  can  be  had  by  constructing  a graph  of  these 
relationships,  and  interpreting  it  as  though  the  cushion  pressure,  cush- 
ion power  and  total  drag  coefficient  were  independent  of  speed. 

This  is  done  in  Figure  A-2,  a logarithmic  graph  of  ideal  specific 
power  versus  the  speed  parameter  V^/^Sp.  The  speed  in  knots  equals 
approximately  17.2  /q"'  for  standard  air  density.  We  consider  V.IJKp  to 
be  nond imens ional , being  merely  a short  way  to  write 


1.69 


Ap 


/0.00238\ 


V 


)} 


On  this  graph,  Equation  [A-l]  gives  the  solid  lines  of  slope  -1;  Equa- 
tion [A-2  ] gives  the  solid  lines  of  slope  +2;  and  Equation  [A-3]  gives 
the  dashed  curves,  for  total  ideal  specific  power,  which  become  asymp- 
totic to  the  cushion  specific  power  lines  at  low  speeds  and  asymptotic 
to  the  propulsion  specific  power  curves  at  high  speeds. 

For  reasons  which  will  become  clear  in  the  succeeding  paragraphs, 
we  have  substituted  an  "ideal  cushion  power  parameter"  ^ (which  is 

identically  equal  to  D for  the  plenum)  in  place  of  the  discharge  coeffi- 

c 

cient  in  Figure  A-2. 
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The  discharge  coefficient  for  a simple  plenum  Lg  typically 

about  0,6  in  magnitude,  the  extreme  range  of  poss Lbilitiea  being 

about  0.5  to  1.0.  Thus,  for  purposes  of  a very  broud  look  at  vehicle 

performance,  we  can  consider  K constant. 

^ > *■ 

Figure  A-2  shows,  then,  the  dominant  influence  on  performance  of 

the  gap  area  ratio  S /S,  and  of  the  drag  coefficient  C . If  we  regard 

8 * 

each  of  the  dashed  curves  as  representing  a specific  design,  then  the 

minima  of  these  curves  are  "optimum"  points,  in  the  sense  of  providing 

minimum  fuel  consumption  during  a given  trip.  Comparing  different 

"designs"  on  this  basis,  the  following  (approximate)  trends  can  be  traced 

on  Figure  A-2: 


1.  With  air  gap  ratio  S /S  fixed,  the  optimum  speed  decreases  with 

8 

- 

increasing  drag  coefficient  in  proportion  to  ' ; and  the  minimum 

1 / 3 

specific  power  increases  as  ' 


2. 


With  drag  coefficient  fixed, 


the  optimum  speed  increases 


with  increasing  air  gap  ratio  in  proportion  to  /s^/S^1/3;  and  the 
minimum  specific  power  increases  in  proportion  to  ^S^/S j 3/3  . 


(These  interpretations  depend  upon  the  cushion  power  and  drag 
coefficient  remaining  nearly  constant  in  the  vicinity  of  the  optimum 
point.  It  will  be  shown  later  that  this  is  usually  approximately  valid, 
provided  V^Z./Sp”  ^ 10,  except,  in  the  case  of  overwater  operation,  near 


the  "hump  speed,"  and  provided  the  cushion  pressure  is  "moderate."  The 
meaning  of  "moderate  cushion  pressure"  will  be  clarified  shortly.) 

Further  discussion  of  Figure  A-2  will  be  postponed  until  we  have 
Introduced  the  peripheral  jet. 


PERIPHERAL  JET  CRAFT 


For  the  plenum  craft,  we  found  the  cushion  power  to  be  proportional 
to  the  product  of  the  3/2-power  of  cushion  pressure  and  the  discharge 
coefficient , 


S 

g 


r 


It  is  possible  to  reduce  the  discharge  coefficient  very  substan- 
tially, by  discharging  the  air  from  the  compressor  at  the  periphery  of 
the  cushion,  instead  of  directly  Into  the  cushion. 


cjo 


— - p 


Plenum 


Peripheral  Jet 


This  reduction  of  the  discharge  coefficient  is  accomplished,  however, 
at  the  expense  of  increasing  the  total  pressure  of  the  compressor  output 
(in  psfg)  from  Ap  to  pt,  pfc  > Ap.  The  cushion  power  thus  Is  proportional 


Pc,l  “ ^ Pt^c 


Rewriting  this  in  the  form 


comparing  with  the  expression  for  the  plenum  given  previously,  and  noting 
further  that  Ap  ■ pf  for  the  plenum,  it  is  found  to  be  convenient  to 


write 


Pc,i  “ ^ Ap  Kc,i 


C,i  Ap/pt 


Plenum 


Peripheral  Jet 
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In  Section  B,  we  shall  consider  in  detail  the  relationships  which 

determine  the  quantities  j?  , Ap/p  j and  their  quotient  K for  the 

c t C j 1 

peripheral  Jet.  For  the  present,  let  us  presuppose  those  results  by 
noting  that,  in  practical  design,  K . is  likely  to  be  approximately 

the  same  (perhaps  20  percent  less)  for  the  peripheral  Jet  ns  for  the 
plenum  craft. 

Other  aspects  of  the  advantages  and  disadvantages  oi  the  periph- 
eral jet  versus  the  plenum  craft  will  also  be  discussed  in  Sections  B 
and  C.  For  purposes  of  our  very  broad  view  of  vehicle  performance, 
however,  it  appears  that  the  two  types  are  about  equivalent.  In 
particular.  Figure  A-2  applies  equally  as  well  to  the  peripheral  jet 
as  to  the  plenum  craft,  as  does,  also,  the  discussion  of  Figure  A-2, 
which  was  carried  out  for  the  plenum  craft.  That  is,  the  minimum 
specific  power  and  optimum  speed  follow  approximately  (with  the  same 
qualifications  pointed  out  previously)  the  proportionality  relation- 
ships : 


C x/a 
D 


cc 

opt 


The  subscript  "opt"  denotes  "optimum"  in  the  aense  of  minimizing  the 
specific  power,  thus  minimizing  the  fuel  consumed  per  mile. 

It  is  very  enlightening  to  trace,  qualitatively,  one  of  the  main 
trends  of  Hovercraft  development  ns  in  the  following  graph  which  should 
be  compared  with  Figure  A-2. 


The  SR.N1  was  allowed  to  have  a relatively  high  drag  coefficient 

C for  reasons  of  econonilcal  construction,  and  had  an  air  gap  area  S 
“1  8 

judged  Just  sufficient  to  avoid  water  contact  under  "normal"  conditions 

This  resulted  in  a rather  large  air  gap  ratio  (S  /S)1  , because  the 

8 

SR. Nl  was  so  small. 

The  SR.N2  was  more  streamlined,  giving  originally  a smaller  drag 

coefficient  C . In  the  original  design,  the  air  gap  height  ("daylight 
Da 

clearance)  h was  about  the  same  as  for  SR.Nl,  but  the  air  gap  area 
ratio  (S  / S)3  was  smaller,  because  the  SR.N2  was  larger.  (For  a given 

O 

planform  shape  and  given  daylight  clearance  h,  the  ratio  S /S  is  in- 

8 

versely  proportional  to  the  linear  dimension  of  the  planfarnu) 


1A. 


However,  before  SR.N2  was  put  In  service,  We3tlnnd  changed  the 
design,  introducing  flexible  trunks  and  reducing  the  air  gap  ratio 
substantially  to  (Sg/S)a'  # xhe  daylight  clearance  was  no  longer  large 

enough  to  avoid  water  contact  under  "normal"  conditions,  so  the  total 
drag  coefficient  under  average  expected  wave  conditions  was  increased 
substantially  ^to  , on  account  of  the  additional  hydrodynamic 

drag. 

The  SR.N4  (projected  150-ton  channel  ferry)  might  have,  according 
to  information  given  in  the  press,  about  the  same  daylight  clearance,  h; 
but  again,  a substantial  decrease  in  air  gap  area  ratio  ^to  (S  /S)4^, 
because  of  the  increased  planfortn  dimensions.  The  average  total  drag 
coefficient  might  go  either  up  or  down,  relative  to  SR.N2,  depending 
on  how  the  respective  average  operating  wave  conditions  were  defined. 

In  the  qualitative  graph  presented  above,  it  has  been  assumed  that 


It  should  be  stressed  again  that  the  above  graph  is  only  quali- 
tative, and  undoubtedly  contains  discrepancies  as  to  the  detailed 
relationships  between  specific  craft.  The  general  trend,  however, 

toward  ever  smaller  gap  area  ratios  S /S,  ever  smaller  specific  power, 

8 

and  ever  smaller  optimum  nondimensional  speeds  is  correctly  represented. 

Also  correctly  represented  is  the  sharp  break  in  the  trend  of 

development  which  occurred  with  the  introduction  of  flexible  trunks. 

One  should  not  consider  that  flexible  trunks  caused  this  break  in  the 

trend;  but  rather,  that  they  permitted  it.  There  was  a substantial 

gain  In  operating  economy  to  be  had  by  reducing  the  air  gap  area  S , 

8 

even  at  the  expense  of  greatly  increased  contact  with  the  surface  (and 
hence  increased  drag)  and  decreased  (nondimensional)  speed.  Flexible 
trunks  made  this  decreased  gap  area  feasible  from  the  standpoint  of 
structural  integrity,  structural  weight,  and  motions  (passenger  comfort). 
(Decreasing  the  air  gap  also  has  other  advantages,  besides  reduced 
specific  power,  as  we  shall  see  in  the  succeeding  sections;  in  partic- 
ular, reduced  size  and  weight  of  machinery  and  ducting.) 


Of  course,  specific  power  Is  not  the  only  performance  parameter 
which  Is  Important.  There  are  three  parameters  which  are  of  paramount 
importance:  specific  power,  ratio  of  empty  weight  to  gross  weight,  and 

speed  (dimensional).  Every  successful  vehicle  Is  the  result  of  a favor- 
able balance  between  these  three  parameters. 

The  unfavorable  trend  toward  lower  nondimens lonal  speeds,  as  the 
specific  power  is  favorably  reduced,  has  tended  to  be  offset  by  increased 
cushion  pressure  so  that  the  dimensional  speed  has  not  decreased.  Actu- 
ally, a gradual  trend  toward  higher  actual  speeds  with  larger  vehicles 
is  justified,  as  we  shall  see  In  later  sections.  In  fact,  If  Figure  A-2 
were  universally  valid,  we  would  not  have  to  worry  very  much  about  speed. 
We  could  accept  almost  any  design  change,  within  reason,  which  improved 
the  specific  power;  and  then  win  back  the  desired  speed  by  increasing  the 
cushion  pressure.  (The  structural  problems  generally  tend  to  become 
more  tractable  also,  with  an  increase  in  cushion  pressure.)  However, 
among  the  various  qualifications  placed  on  Figure  A-2  was  the  requirement 
that  the  cushion  pressure  be  "moderate."  It  is  now  necessary  to  clarify 
this  requirement,  which  has  to  do  with  the  wavemaking  drag. 

The  practical  validity  of  Figure  A-2,  and  of  the  various  arguments 
we  have  carried  out  in  connection  with  this  figure,  depends  upon  the 
following  conditions:  In  the  vicinity  of  the  optimum  point, 

a.  The  total  drag  is  varying  nearly  in  direct  proportion  to  V^s. 

b.  The  cushion  power  is  nearly  constant. 

Since  some  of  the  components  of  the  total  drag  are  definitely  following 
decidedly  different  laws  from  the  V^3  law,  Figure  A-2  is  meaningful  ^mly  . 
if  these  components  are  a relatively  small  part  of  the  total  drag,  in 
the  vicinity  of  the  optimum  speed. 

The  wavemaking  drag,  in  particular,  varies  in  proportion  to  a 

negative  power  of  at  speeds  above  "hump  speed";  therefore.  Figure 

A-2  does  not  apply  at  all  if  the  wavemaking  drag  is  very  significant,  in 

the  vicinity  of  the  optimum  speed.  Moreover,  for  a given  planform  shape 

and  size,  the  wavemaking  drag  D divided  by  the  gross  weight  W varies 

w 

directly  with  the  cushion  pressure  Ap: 


Therefore,  when  we  say  Figure  A-2  is  restricted  to  "moderate"  cushion 
pressures,  we  mean  specifically  that  it  is  restricted  to  cushion  pres- 
sure sufficiently  moderate  that  the  wavemaking  drag  plays  only  a minor 
role  in  the  determination  of  the  optimum  performance  point. 

The  Westland  family  of  Hovercraft  SR.Nl  through  SR.N4,  which  were 
discussed  above,  and  also  SR.N5  and  SR.N6,  which  are  reversions  to  a 
smaller  size,  satisfy  the  condition  of  "moderate"  cushion  pressure 
reasonably  well.  However,  it  would  seem  that,  projecting  that  line  of 
development  into  the  future,  to  larger  and  larger  vehicles,  the  wave- 
making drag  would  become  more  and  more  important,  and  optimum  performance 
would  be  determined  less  and  less  adequately  by  the  simple  cushion  power 
versus  propulsion  power  tradeoff  represented  in  Figure  A-2. 

Rather  than  attempt  to  follow  the  subtle  changes  likely  to  take 
place  in  that  line  of  development,  however,  it  will  be  simpler  for  our 
present  purpose  to  turn  to  a different  line  of  development  in  which  we 
already  find  the  cushion  power  playing  a minor  role,  while  the  wavemaking 
drag  plays  a very  major  role. 

CAPTURED  AIR  BUBBLE  (CAB),  OR  SIDEWALL  GEM 

Many  of  the  early  ideas  and  inventions  of  Ground  Effect  Machines 
included  the  idea  of  minimizing  the  air  gap  by  use  of  slender  sidewalls, 
which  extended  into  the  water  to  block  the  escape  of  air  from  the 
cushion. 


Plenum-Type  Sidewall  Ground  Effect  Machine 

In  the  intensive  early  exploration  of  GEM's  (from  1957  to  1961, 
say),  the  sidewall  GEM  received  its  due  share  of  attention;  but  most 
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Investigators  quickly  concluded  (perhaps  prematurely)  that,  because  of 
the  hydrodynamic  drag  of  the  sidewalls,  these  craft  would  be  suitable 
only  for  relatively  low  speeds,  perhaps  of  the  order  of  50  knots. 

Two  separate  trends  have  developed  to  change  this  outlook.  First, 
as  we  have  just  noted,  the  foremost  developers  of  "ordinary"  OEM's  (we 
will  call  them  "full-peripheral"  GEM's  to  denote  more  or  less  uniform 
air  gap  h all  around)  are  finding  it  profitable  to  accept  much  more 
severe  hydrodynamic  drag,  in  the  interest  of  reducing  the  air  gap  area 

S , than  would  have  been  guessed  by  most  experts  a few  years  ago. 

8 

Secondly,  more  recent  research  on  sidewall  GEM's  (especially  under  the 
Captured  Air  Bubble  research  program  at  the  U.  S.  Naval  Air  Development 
Center)  has  produced  greatly  improved  understanding  of  the  sidewall  GEM. 
The  sidewall  GEM,  or  CAB  now  appears  to  be  a very  strong  contender  for 
future  maritime  GEM  applications. 

For  the  full-peripheral  GEM  (at  its  current  stage  of  development, 
with  "moderate"  cushion  pressures),  the  most  essential  thing  to  be  under- 
stood is  the  tradeoff  between  cushion  power  and  propulsion  power.  For 
the  CAB,  the  most  essential  thing  to  be  understood  is  the  tradeoff  between 
the  wavemaking  propulsion  power  and  the  remainder  of  the  propulsion  power. 

Previously  we  dealt  with  the  ideal  total  specific  propulsion  power 


° * c 
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We  must  now  split  this  into  two  parts,  "wavemaking"  and  "other" 
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For  a given  planform  shape,  the  wavemaking  drag  can  be  represented 


W 


functionally;  that  is. 


where  V^JJT  is  the  "speed-length"  ratio  (a  nond linens ional  speed  pa- 
rameter, in  the  sense  that  it  is  derived  from  the  nond iraena ional  Froude 
number  V q/ by  inserting  a constant  value  for  the  acceleration  of 

gravity  g) . (Similarly,  Ap/X  is  nond lmens ional  in  the  sense  of  being 
derived  from  Ap/(\^  X)  by  inserting  a constant  value  for  the  weight 

density  of  water  v ) • 

* • 

For  simplicity  we  shall  consider  only  a single  planform  shape 

(rectangular,  X/b  » 2)  and  the  corresponding  single  graph  of  —■  -r  4^ 

w l 

versus  V^/ JT  (Figure  A-3)  for  purposes  of  this  introductory  discussion. 

When  we  seek  to  represent  the  CAB  on  a graph  like  Figure  A-2,  we 
find  the  problem  slightly  complicated  by  the  addition  of  the  two  new 
variables  Ap/X  and  . However,  one  of  these  is  eliminated  by  re- 

lating the  speed-length  ratio  V^/^/Ap  as  follows: 


- (yv®)  JSpi 


Hence,  only  the  single  new  variable  Ap/X  need  appear  on  the  graph. 
Such  a graph  is  presented  in  Figure  A-4,  which  shows,  for  simplicity. 


only  single  representative  values  of  the  variables  C_  and  K . [ _g)  . 

Do  ’ \ S / 

Figure  A-4  displays  straight  lines  of  slopes  -1  and  +2,  just  like 
Figure  A-2;  but  now,  the  intersection  between  these  lines  is  faired 
by  a more  complicated  curve  of  total  ideal  specific  power,  strongly  in- 
fluenced by  the  wavemaking  drag  "hump."  In  fact,  we  now  find  two  minima 
in  the  "total"  curve  when  the  cushion  is  heavily  loaded.  The  minimum  at 

the  higher  speed  is  the  one  which  we  choose  to  call  the  optimum  perform- 
ance point.  (In  practice,  the  air  gap  ratio  S /S  might  be  reduced  as 

8 

the  speed  is  reduced  to  such  an  extent  that  the  lower  minimum  disappears. 
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and  the  specific  power  Improves  continuously  with  speed  reduction,  at 
speeds  below  "hump  speed,"  In  fact,  the  cushion  power  has  usually  been 
neglected  entirely  from  the  optimization  studies  for  CAB's,  except  that 
a 5 or  10  percent  power  margin  is  added  at  . 'nd.) 

One  sees  In  the  example  of  Figure  A-4  that  an  Increase  In  cushion 
loading  parameter  Ap//  has  a very  slightly  unfavorable  effect  on  the 
value  of  the  minimum  specific  power,  at  the  higher  speed  minimum.  On 
the  other  hand,  an  Increase  in  Ap //  has  a noticeably  favorable  effect 
on  the  (dimensional)  optimum  speed.  However,  at  very  high  loadings,  the 
severe  wavemaking  drag  hump  complicates  the  design  of  the  propulsion  sys- 
tem. 


K 

c. 
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Keep  in  mind  that  Figure  A-4  gives  only  one  example  of  C 


D 


and 


o 

It  is  easy  to  see  that  the  picture  can  change  considerably 


with  changes  In  these  variables.  . 

The  alert  reader  will  now  ask,  Wherein  has  our  preliminary  treatment 
of  CAB's,  as  reflected  In  Figure  A-4,  and  our  preliminary  treatment  of 
full-peripheral  GEM's,  as  reflected  in  Figure  A-2,  revealed  any  funda- 
mental differences  between  the  two  types  of  craft?  We  have  no  answer, 
except  to  admit  that  we  have  found  no  fundamental  differences  at  all! 

In  treating  the  CAB,  we  have  been  compelled  to  bring  in  the  concept  of 
wavemaking  drag,  because 


a.  We  expect  such  small  values  of  cushion  specific  power  that 
the  component  of  propulsion  specific  power  associated  with  wavemaking 
drag  becomes  significant  by  comparison,  and 

b.  We  expect  somewhat  higher  values  of  the  parameter  Ap //  for*  the 
CAB,  which  makes  the  wavemaking  drag  larger  in  absolute  magnitude. 


However,  these  are  differences  of  degree,  not  differences  of  kind.  The 
fact  Is  that  we  might  just  as  well  have  introduced  wavemaking  drag  into 
Figure  A-2  (had  we  not  been  striving  for  simplicity);  and  the  fact  is 
that  Figure  A-4  1s  every  bit  as  valid  for  a full-peripheral  GEM  (with 

K , ( _ S ) B 0.002,  C_  - 0.5  and  rectangular  cushion,  (//b  = 2)  as  it  is 
C»1\S  / °0 
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for  a GAB.  Of  course,  such  a small  value  for  the  air  gap  area  ratio  Is 
not  at  all  representative  of  current  full-peripheral  GEM's,  but  the 
trend  of  development  seems  to  point  somewhat  In  that  direction. 

Note  that,  since  the  cushion  power  Is  such  a small  fraction  of  the 
CAB's  total  power,  there  Is  little  or  nothing  to  be  gained  from  sophis- 
tication In  the  method  of  delivering  air  to  the  cushion.  Designers  will 
probably  choose  to  let  the  air  be  delivered  directly  from  the  compressors 
to  the  cushion,  as  In  the  full-peripheral  plenum  craft  and  as  indicated 
schematically  in  the  sketch  at  the  beginning  of  our  discussion  of  CAB's, 

Also,  referring  to  that  sketch.  It  Is  apparent  that  the  plenum 
walls  at  the  bow  and  stern  will  preferably  be  compliant  with  the  uneven 
surface  (waves)  over  which  the  craft  travels.  Two  of  the  proposed  solu- 
tions are: 

a.  To  spring-suspend  the  bow  and  stem  plenum  walls  mechanically, 
in  such  a way  that  they  will  plane  on  the  wave  surfaces,  as  skiis,  at  a 
favorable  angle  of  attack  for  minimum  drag;  or 

b.  To  use  flexible-fabric  bow  and  stem  plenum  walls  arranged  to 
yield  freely  to  wave  impacts;  In  other  words,  to  fit  the  CAB  with  flexible 
trunks, 

TRENDS 

The  influence  of  flexible  trunks,  and/or  of  sidewalls  and  flexible 
bow  and  stern  seals,  can  be  viewed  at  this  point  as  aspects  of  a con- 
sistent general  trend:  Reduction  of  the  air  gap  area  ratio,  reduction  of 

the  specific  power,  and  reduction  of  the  nondimens ional  speed.  Progress 
to  larger  and  larger  vehicles  can  be  seen  to  result  In  a further  continu- 
ation of  this  same  trend. 

This  trend  is  indicated  in  Figure  A-5,  which  Identifies  the  domain 
of  current  GEM's  to  be  specific  powers  of  the  order  of  0.2  to  0.3  (roughly 
the  same  as  helicopters)  with  speed-pressure  ratios  of  the  order  of  5 to 
10;  and  implies  a future  extension  of  this  domain  into  much  lower  specific 
powers  and  somewhat  lower  nondimens lonal  speeds.  It  Is  emphasized  again 
that  reduction  of  nondimens ional  speed  does  not  necessarily  Imply  lower 
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actual  speed  In  knots.  In  fact,  to  the  extent  that  the  reduction  in 
nondimenaional  speed  results  from  progress  to  larger  vehicles,  the  speed 
in  knots  will  be  found  to  increase. 

A vertical  dashed  line  on  Figure  A-5  marks  the  value  (17.2)  of 
the  speed  parameter  at  which  the  air  dynamic  pressure  associated  with 
forward  speed  equals  the  quantity  W/S  (cushion  loading,  wing  loading,  or 
disc  loading  for  the  various  vehicles  represented)  at  sea  level  standard 
conditions.  To  the  right  of  this  line  lie  ground  effect  vehicles,  such 
as  the  "Ram  Wing"  and  "Wing  in  Ground  Effect"  (WIG),  which  are  sustained 
by  pressures  induced  by  the  forward  motion  rather  than  by  a compressor- 
fed  cushion.  It  was  once  thought  that  the  future  evolution  of  GEM's 
might  proceed  in  this  direction;  that  is,  toward  more  aircraft-like 
vehicles.  However,  the  current  trend  appears  to  be  definitely  toward 
more  ship-.like  GEM's,  as  indicated  in  Figure  A-5. 


Aerodynamics  Laboratory 
David  Taylor  Model  Basin 
Washington,  D.  C. 

March  1966 
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Figure  A-l  - Fuel  Consumed  in  Various  Distances  as  Function  of  Total  Specific  Power 


Ideal  Specific  Power 


Definitions  and  Expected  Ranges 

Plenum: 

Kc,t  ■ . °-5  < Kcl  < 1.0 

Peripheral  Jet: 

K , s & * A P/P.  , 0.4  < K . < 1.0 

x c c C,  1 

[W  = Ap  S ; pa  = 0.00238  s lugs/ ft3  3 


Figure  A- 2 - Ideal  Specific  Power  Versus  Speed  Parameter  as 
Function  of  Gap  Area  Ratio  and  Total  Drag  Coefficient 
(Effect  of  forward  speed  on  cushion  power  is  neglected.) 


Speed-Length  Ratio,  V^/^/T 


Figure  A- 3 - Wavemaking  Drag  for  a Rectangular  Cushion  of 
Length/Beam  Ratio  l/b  - 2 

(According  to  Newman  and  Poole,  DTMB  Report  1619,  March  1962) 
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Foreword 

This  report  is  based  on  a lecture  series  presented  by  the 
authors  at  the  von  Karman  Institute  for  Fluid  Dynamics,  Rhode- 
Saint-Genese,  Belgium,  in  May  1965;  and  at  the  University  of 
Maryland,  College  Park,  Maryland,  in  July  1965.  The  lectures  were 
prepared  under  the  joint  auspices  of  the  David  Taylor  Model  Basin 
and  the  Naval  Air  Development  Center.  They  were  presented  in 
Belgium  under  the  joint  sponsorship  of  the  von  Karman  Institute 
and  the  Advisory  Group  for  Aerospace  Research  and  Development 
(AGARD);  and  in  Maryland  under  the  sponsorship  of  the  Assistant 
Secretary  of  the  Navy  for  Research  and  Development. 

The  revised  lectures  will  be  presented  as  follows: 

A.  Introductory  Survey 

B.  Air  Cushion  Mechanics 

C.  Internal  Aerodynamics 

D.  Drag 

_g-« Opr  jmi  T’at-i  on— 1 -^FM 

F.  Cushion  Contributions  to  Stability 

G.  Seakeeping 

H.  Performance  Summary 


NOTATION 

cushion  length,  ft 
cushion  beam,  ft 

cushion  area,  ft2,  measured,  in  plan  view,  to  outer 
edge  of  nozzle  exit 

cushion  pressure,  psfg 

total  pressure  of  air  supply  to  cushion  at  nozzle  exit 
or  cushion  entry,  psfg 

flight  velocity  in  earth  axes,  ft/sec 

flight  velocity  in  earth  axes,  knots 

air  density,  slugs/ft3 

water  density,  slugs/ft3 

reference  air  dynamic  pressure  Vq2/2  j , lb/ft2 

reference  water  dynamic  pressure  ( p^  Vq2/2^  , lb/ ft2 

local  air  velocity,  ft/sec 

local  dynamic  pressure  (pu2/2),  lb/ft2 

daylight  gap  area  (cushion  perimeter  times  mean  daylight 
gap  (hC), ft2) 


"daylight"  gap,  feet  (local  value  or  mean  value  as 
appropriate) 

cushion  perimeter,  feet 

peripheral  jet  efflux  angle,  measured  positive  inward 
from  vert ical, degrees 

nozzle  thickness,  feet 

nozzle  thickness  parameter ^ ^ (1  + 

total  air  volume  flow  rate  through 

! 

cusk  discharge  coefficient  / _~—’rr  | 


sin  6)^ 

cushion  system,  ft3/s 

Q \ 


jet  reaction  coefficient 

gross  weight,  pounds 
cushion  loading  ^ j > lb/ft2 

external  aerodynamic  lift,  pounds 

external  aerodynamic  drag,  pounds 

external  aerodynamic  lift  coefficient  ^Lg/qaSj 

external  aerodynamic  drag  coefficient  ( D /q  S\ 

\ e a / 


total  shaft  power  supplied  to  cushion  system,  lb-ft/sec 

total  shaft  power  supplied  to  propulsion  system,  lb-ft/sec 

total  shaft  power  Pc  + j , lb-ft/sec 

ram  recovery  efficiency 

duct  efficiency 

compressor  efficiency 

internal  efficiency 


iv 
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SUMMARY 

The  energy  and  momentum  relationships  governing  air  cushion 
performance  are  reviewed.  The  exponential-theory  equations  are  rec- 
ommended for  calculation  of  cushion  pressure  and  jet  reaction,  and  a 
modified  equation  Is  proposed  for  calculation  of  the  volume  flow  rate. 

The  question  of  off-design  cushion  performance  is  discussed.  It 
is  pointed  out  that  off-design  performance  is  of  paramount  practical 
importance  and  that  the  widely  held  notion  that  peripheral  jet  cushions 
are  more  efficient  than  simple  plenum  cushions  is  not  necessarily  vaild 
in  this  context. 

The  danger  of  drawing  erroneous  conclusions  from  the  cushion 
performance  equations,  which  (except  for  the  plenum)  do  not  apply  to 
off-design  operation,  is  emphasized. 

INTRODUCTION 

There  have  been  almost  as  many  formulations  of  the  fluid  dynamic 
relationships  governing  the  behavior  of  air  cushions  as  there  are  fluid 
dynamicists  working  in  this  field.  The  air  cushion  poses  intriguing  and 
important  problems  which  — like  most  problems  of  fluid  mechanics  — can 
become  extremely  complex  if  one  seeks  to  formulate  them  exactly. 

The  fact  is,  however,  that  experimental  measurements  of  the  im- 
portant quantities  characterizing  air  cushion  performance  fall  into  a 
rather  simple  pattern,  compatible  with  comparatively  simple  theoretical 
treatments . 

We  will  consider  three  separate  formulations:  The  "Thin-Jet  Theory 

which  is  of  limited  practical  value,  but  provides  an  easily  understood 
introduction  to  some  main  ideas  of  the  peripheral  jet;  the  "Exponential 
Theory,"  which  comes  closest  of  all  the  simple  theories  to  providing 
solutions  fully  satisfactory  for  design  application  (but  fails  to  cor- 
rectly predict  one  of  the  essential  quantities  in  the  case  of  very  thick 
jets);  and  Plenum  Theory.  Finally,  we  will  use  the  results  of  the 
Plenum  Theory  to  arrive  at  a modification  of  the  Exponential  Theory  re- 
sults, providing  simple  formulas  which  are  suitable  for  practical  design 
application. 


#PREC£Dl!75 


LAM5C.N0T  F 


All  theoretical  considerations  will  be  pursued  under  the  as- 
sumption that  air  will  behave  as  an  Incompressible,  Invlscid  fluid. 
Moreover,  we  will  employ  throughout  the  approximation  that  the  airflow 
at  the  cushion  periphery  Is  essentially  two-dimensional. 

THEORETICAL  FORMULATIONS 

THIN- JET  THEORY 

Consider  an  air  cushion  bounded  top  and  bottom  by  impervious 
surfaces , and  separated  from  the  outside  atmosphere  by  a peripheral 
jet  ("air  curtain"),  as  shown. 


Neglecting  mixing  between  the  jet  and  Its  surroundings,  and 
assuming  the  cushion  to  be  in  static  equilibrium  at  pressure  Ap  (no 
air  entering  or  leaving  cushion),  It  will  be  seen  that  the  jet  must 
curve  and  become  tangent  to  the  ground,  as  shown  in  the  sketch. 

Bernoulllrs  equation  applied  to  the  jet  is 

p + q “ Pt 


2 


where  pt  is  the  total  pressure  (assumed  to  be  uniform  throughout  the 

Jet),  p is  local  static  pressure,  and  q is  local  dynamic  pressure  (pu2/2 , 
where  u is  the  magnitude  of  the  local  velocity).  (All  pressures  are 
"gage,"  measured  relative  to  the  outside  atmosphere.) 

Thus,  the  dynamic  pressure  along  the  outer  free  boundary  of  the 
jet  (p  = 0)  is  q = p^;  and  along  the  inner  free  boundary  (p  = Ap),  it 
is  q = pt  - Ap. 

Now  if  we  introduce  the  assumption  * • 

Ap  « Pt 

it  will  be  seen  that 


q = 


Pt 


throughout  the  jet;  and  the  magnitude  of  the  jet  momentum  per  unit 
cushion  periphery  at  nozzle  exit 


= J pu2  dy  = 2 t = 2 P(;  t 

o 

will  also  be  nearly  uniform  along  the  free  path  of  the  jet. 

The  horizontal  pressure  force  Ap*h  per  unit  cushion  periphery  must 
be  balanced  by  the  change  in  horizontal  component  of  jet  momentum: 


or 


where 


Ap *h  = jL  sin  0 + j2  = 2 pt  t (1  + sin  0) 


x « 1 


X 


t 

h 


(1  + sin  0) 


Lb-i3 


The  nondimens ional  quantity  x is  called  the  "nozzle  thickness  param~ 
eter."  We  are  clearly  justified  in  replacing  the  original  assumption 
(Ap  « Pt)  by  the  condition  x « 1. 


3- 


Equal  in  importance  to  the  cushion  pressure  ratio  Ap/p,.  is  the 
discharge  coefficient  .0,  which  relates  the  volume  rate  of  flow  Q, 
through  the  cushion  system,  to  the  cushion  pressure.  We  define 


A s 


c sgyr-5 


where  S (which  equals  hC  in  che  case  under  consideration)  is  the  area 

s 

of  the  "daylight"  gap  at  the  edge  of  the  cushion. 

We  have  already  seen  that,  if  x « 1,  the  jet_  velocity  at  the  nor.- 


uniform 

and 

of 

magnitude 

Q = 7 

1 P 
P 

> i 

1 

> 

J3  = / 

'it 

t 

c 

Ap 

h 

Equation 

[B- 

i], 

we 

obtain 

1 

f! 

c 1 + 

s in 

"e 

■v  2 

> 

« 1 


« 1 


[B-2] 


We  have  already  found  an  expression  for  the  jet  momentum,  which 
immediately  provides  an  approximate  solution  for  the  jet  reaction 
coefficient 


C.  = 
J 


Ap  t 


. 1 


x 


x « l 


Lb-33 
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EXPONENTIAL  THEORY 


For  values  of  the  nozzle  thickness  parameter  x of  Interest  for 
practical  design  (say,  x > 0.2),  the  thin-jet  theory  does  not  give 
satisfactory  results.  In  order  to  obtain  practical  design  formulas, 
it  is  necessary  to  take  account  of  the  variation  in  static  pressure 
across  the  jet.  Different  investigators  have  done  this  in  many 
different  ways,  including  direct  postulation  of  the  static  pressure 
distribution,  or  postulation  of  the  radius  of  curvature  R(y)  of  the 
jet  streamlines  at  the  nozzle  exit,  whence  the  pressure  variation  ear. 
be  calculated  from  the  relation 


iE  = _ jjuf.  = _ - S 

dy  R R 

The  most  successful  of  the  simple  theories  of  this  kind  was  due 
to  Stanton-Jones  and  Els  ley,  who  assumed  that  an  adequate  approximation 
could  be  obtained  by  setting 


One  can  then  write,  from  the  Bernoulli  equation  and  the  above  relations: 


|e  = - i - 3.  (1  + sln  Q) 
dy  dy  h 


where  x = ^ (1  + sin  0),  as  before.  This  equation  is  readily  integrated 
to  obtain 

log  q = 2x  ^ + Constant 
e t 


5 


where  the  constant  can  be  determined  from  the  boundary  conditions 


yielding 


q (0)  = Pt  - Ap 


q (t)  = pt 


q = (Pt 


y 

T 
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Continuing,  the  discharge  coefficient  can  be  calculated 


c 4 pAph 


Substituting  for  ^ from  Equation  [B-4]  gives,  after  some  reduction 


U3  = q t / tanh  ^ 

c 1 + sin  0 Ay/  2 


[B-5] 


(Note:  Do  not  use  this  equation  if  x > 0.4.  See  Equations 

[B-7]  and  [B-5 1 ] 


The  jet  reaction  coefficient  is  given  by 


cj  ■ ^ jV* + >>>dy  - ^ J <pt  + <•> dy 


P*  1 ct| 


2x^1  . 
e t dy 


I 


Note  that  Equations  [B-4 ] through  [B-6]  agree  exactly  with  Equations 
[B-l]  through  [B— 3 ] in  the  limit  x -*  0. 

It  is  found  that  Equation  [B-4  ] agrees  remarkably  well  with  experi- 
ment in  the  whole  range  x > 0.2,  provided  adequate  account  is  taken  of 
non-uniformities  in  total  pressure  over  the  nozzle  exit.  (In  the  range 
x < 0.2,  the  equation  is  valid  for  the  ideal  fluid  considered;  but,  for 
real  fluids,  mixing  effects  become  important  in  this  range.) 

Equation  [B-5],  on  the  other  hand,  is  found  to, agree  with  experiment 
in  a more  limited  range,  0.2  < x < 0.4,  say.  The  Plenum  Theory  will 
provide  the  basis  for  a less  restricted  formula  for 

Experimental  data  available  from  which  to  draw  conclusions  with 
regard  to  the  validity  of  Equation  [B-6]  are  limited.  The  equation  is 
almost  certainly  reliable  in  the  range  x < 0.4,  and  it  seems  quite  likely 
to  be  valid  for  most  practical  purposes  in  the  whole  range  0 < x < <*>. 

In  any  event,  it  will  be  shown  later  that  the  jet  reaction  coefficient 
plays  a rather  small  role  in  most  practicaal  design  problems. 

PLENUM  THEORY 

It  is  clear,  on  physical  grounds,  that  the  peripheral  jet  craft  and 
the  plenum  craft  become  indistinguishable,  from  the  standpoint  of  fluid 
mechanics,  when  the  former's  nozzle  thickness  parameter  x becomes  very 
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large.  The  plenum  craft  could  be  defined  as  the  limiting  form  of  the 
peripheral  jet  craft,  in  the  limit  x-*00  • 


p *=  Ap 


p « 0 


^c  h 


/ / / / / / / 


Equation  [B-4 ] gives 


(£)-  1-« 


and  Equation  [B-6]  gives 


Lira  C.  = 1.0 

X-.00  J 


The  results  in  both  cases  agree  with  our  intuition  and  experience. 
However,  Equation  [B-5]  gives 


Llm  jy , 


1 + sin  9 


which  does  not  agree  with  experience,  except  in  the  isolated  case 
0 = 90°.  It  is  therefore  to  the  discharge  coefficient  .0  that  we  will 
direct  our  attention. 

Unfortunately,  the  seemingly  simple  fluid  dynamic  problem  of  dis- 
charge through  a sharp-edged  orifice  presents  rather  formidable  mathe- 
matical difficulties  when  one  attempts  direct  solution. 


I 

I 


Classical  hydrodynamics  provides  solutions  in  the  following 
discrete  cases  t 

Discharge  Coefficient  for  Plenum  Craft 


6,  deg 

90 

0 

-90 

c 

1 

TT 

1.0 

2 

TT  + 2 

Further  solutions  are  not  readily  available.  F.  Ehrich  has  presented 
a numerical  method  of  calculation  (Journal  of  the  Aerospace  Sciences, 
November  1961),  but  numerical  results  of  use  in  the  present  problem  were 
not  included.  It  is  understood  that  A.  Gabbay  has  obtained  a general 
solution  (University  of  London  Thesis,  1960),  but  this  is  not  available 
to  the  present  author  at  this  time. 

We  will  therefore  proceed  to  select,  rather  arbitrarily,  a formula 

which  is  consistent  with  the  facts  at  hand.  In  addition  to  the  discharge 

coefficients  presented  above,  some  additional  inferences  can  be  drawn  from 

consideration  of  the  pressure  forces  acting  on  the  orifice  lip.  These 

considerations  are  discussed  in  Section  F.  We  will  presuppose  them  here 

by  stating  that,  if  the  formula  for  & is  written  iif  the  form 

c 


where  f ( 0)  is  a function  yet  to  be  determined,  it  can  be  argued  on 
physical  grounds  that 
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which  satisfies  these  conditions,  as  well  as  the  facts  presented  in 
the  table,  is 


f(9)  “ " +-l  (1  + sin  0)  - sin  9 cos  9 
rr  - l 

We  will  therefore  adopt  the  relationship 


[B-7] 


There  is  no  question  of  the  adequacy  of  Equation  [B-7  ] for  practi- 
cal design  purposes.  The  empirical  nature  of  this  equation  should  be 
borne  in  mind,  however,  if  one  is  tempted  to  apply  it  to  analysis  of  more 
subtle  academic  questions.  For  example,  it  will  be  shown  in  Section  F 
that  the  quantity  2 - 1 has  a certain  significance.  While  we  can  be 
confident  that  the  percentage  error  incurred  in  predicting  £ c itself  by 
Equation  [B-7]  is  small,  a substantial  percentage  error  might  be  incurred 
in  predicting  - 1.  (Care  has  been  taken,  however,  to  avoid  any 
order-of-magnitude  error  in  this  particular  quantity.) 


MODIFIED  EXPONENTIAL  THEORY 

The  Exponential  Theory  gives  results  (Equations  [B-4  ] through  [B-6]) 
which  are  suitable  for  design  application,  except  that  Equation  [B-5  ] 
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determined  from  the  Plenum  Theory  that  the  correct  limits  should  be  ap- 
proximately as  given  by  Equation  [B-7  ].  It  will  be  convenient  to  have 
a new  expression  for  discharge  coefficient  which  agrees  with  Equation 
[B-5]  for  small  and  moderate  values  of  x,  but  approaches  agreement  with 
Equation  [b— 7 ] as  x — Such  an  expression  (Equation  [B-51  ])  is  given 
below  as  part  of  a complete  set  of  air  cushion  relationships  suitable 
for  design  application: 

If  x > 0.2: 


J3  = J3  / tanh 
c c / 


[-V1 


where 


jy  = ft 

c x-»®  c 


i \l  + v 

L rr 


cos  0 

> 

■(I  + sin  0)  - sin  0 cos  0. 


[B-7] 


-2x  + 2x 


[B-6] 


These  relationships,  which  will  provide  the  basis  for  most  of 
our  subsequent  considerations  of  the  air  cushion,  are  presented  graphi- 


cally in  Figures  B-l  through  B-3,  along  with  the  quantity 


Ap/pt 


which  has  a special  significance  in  connection  with  the  cushion 
power  — our  next  topic  of  discussion. 

Note  particularly  the  fact  that  we  no  longer  need  to  draw  any 
sharp  distinction  between  cushions  of  the  plenum  type  and  of  the  pe- 
ripheral jet  type.  The  above  relationships  yield  cushion  properties 
which  gradually  approach  those  of  the  classical  plenum  type,  as  x is 
made  larger  and  larger. 


PERFORMANCE  CHARACTERISTICS 
TOTAL  LIFT  AND  CUSHION  POWER 

The  total  lift  must  equal  the  gross  weight  W,  In  an  equilibrium 
condition.  It  Is  composed  of  the  sum  of  the  base-pressure  reaction, 
the  jet  reaction,  and  the  external  aerodynamic  lift. 

W = Ap  S + (Cj  - 1)  Ap  SN  cos  9 + CL  qfl  S 
or 

W = Ap  s[l+  (Cj  - 1)  cos  0 + CL  ^ J [B-8 ] 

where  SM  is  total  nozzle  area  and  q is  free-stream  dynamic  pressure 
of  the  ambient  air. 

The  cushion  power  P^  is  the  shaft  power  supplied  to  the  cushion 
system  compressor (s) . In  the  absence  of  any  losses,  this  would  simply 
be  the  rate  of  work  required  to  raise  the  total  pressure  of  air  at 
volume  flow  rate  Q from  qa  (free-stream  total  pressure)  to  pfc  (nozzle- 
exit  total  pressure);  that  is, 

Pc  = Q (pt  * qa> 

However,  there  will  inevitably  be  losses  in  the  internal  system.  These 
will  be  discussed  at  some  length  in  Section  C.  For  the  present,  we 
will  simply  introduce  and  apply  the  following  efficiency  parameters,  to 
be  defined  in  Section  C: 

Recovery  efficiency, 

Duct  efficiency, 

Compressor  efficiency,  7^. 

Internal  efficiency,  1]^  = 7]^ 
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The  cushion  power  becomes : 


P 

c 


Recalling  the  previous  definition 


c 


S 

g 


we  have,  after  some  reduction 


[T 


P = Ap  / - Ap  S 
c V p g 


C . _ a V CL 

Ap/pt  ' 'int  c T)c  Ap 


[B-9] 


^ c 

The  significance  of  the  quantity  , shown  graphically  in  pig- 


ure  B-3,  is  now  apparent.  For  design  conditions  common  to  most  of  the 
important  current  vehicles;  namely, 


Is 

S 


« 1 


1 


the  cushion  power  is  almost  directly  proportional  to  this  quantity. 

If,  subject  to  these  conditions,  we  were  to  define  an  "optimum  jet 

geometry"  providing  minimum  cushion  power  for  given  weight  W and  cushion 

area  S,  Figure  B-3  would  show  this  optimum  to  correspond  roughly  to  the 

geometry:  9 90°,  x = 0.7.  — = 0.35.  However,  a rather  large  number 

h 

of  practical  considerations  intervene  to  render  any  such  simple  concept 
of  an  "optimum  geometry"  practically  meaningless. 


Equation  [B-9]  has  the  weakness  that  it  contains  the  quantity 
Ap,  the  cushion  pressure,  which  is  not  really  a primary  design  variable, 
although,  as  we  have  already  noted,  Ap  is  usually  approximately  equal 
to  the  primary  variable  "cushion  loading"  w,  defined  as  follows: 

_ W 

W = — 

s 

We  can  rectify  the  situation  by  substituting  for  Ap  from  Equation  [b-3 j, 
but,  in  so  doing,  we  will  introduce  the  important  §impi  .-'/ing  approxi- 


mation 


iP  s „(i  - cL^)  , / « i 


tB-8'  ] 


(The  justification  for  the  condition-^  « 1 is  as  follows:  Referring 

to  Equation  [B-8],  it  is  evident  that  the  approximation  [b-81]  is  valid 


provided 


Noting  that 


(C,  - 1)  cos  0 « 1 
J *> 


SN  5 tC 


S = hC 
g 


1 + sin  0 


we  see  that  Equation  [B-8' ] is  valid  provided 


.(C  _ n _£  -co?-J  - « i 

(Cj  ' S 1 + sin  0 
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Now  referring  to  Equation  I.B-6],  we  find 


The  ambitious  reader  can  convince  himself  that  the  quantity 

/_  cos  e 

x , C . - 1 ) ■: — ; : ; 

^ j / 1 + sin  0 

is  of  order  unity  or  lower,  under  all  conditions  of  practical  interest 

S 

hence,  the  condition  « 1 is  sufficient.) 


Combining  Equations  [B~9j  and  [B-8  j,  we  obtain 


It  is  convenient  to  rewrite  this  relationship  in  terms  of  the 
nondimensional  "cushion  specific  power": 


3 


where 


S 

« l [b-io] 
s 
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An  easier  insight  into  the  relationship  given  by  Equation  [B-10] 
is  obtained  by  substituting  q = (V,/17.2)s  for  standard  air  density, 

8 K 

and  introducing  a further  approximation,  valid  when  (q  /w)  « 1;  that 

3 

is, 


Pc  _.  17 . 2/w/  l 

wv  s c v,  /yj 

o k v 


pq aM  « i 

|^(Sg/S  « l 


[B-10'  ] 


where  w is  expressed  in  pounds  per  square  foot  and  in  knots.  The 
condition  q /w  « 1 is  not  quite  satisfied  throughout  the  whole 

SL 

operating  range  of  practical  vehicles,  so  Equation  [b-101 ] is  not 
suitable  for  detailed  design;  but  it  correctly  expresses  the  dominant 
behavior  of  the  cushion  specific  power.  The  essence  of  the  truth  is 
expressed  in^ the  not-quite-true  statement,  "The  cushion  specific  power 
is  directly  proportional  to  the  ratio  of  air-gap  area  to  cushion  area, 
and  inversely  proportional  to  the  speed-length  ratio  V^//T." 

The  total  specific  power  is  the  sum  of  cushion  specific  power 
and  propulsion  specific  power 


o o o 


and  is  directly  proportional  to  the  fractional  weight  of  fuel  con- 
sumed per  mile.  (The  factor  of  proportionality  for  ordinary  engines 

p 

is  approximately  1/650.  Roughly  speaking,  then,  if  — = 0.1,  about 

WVq 

10  percent  of  the  gross  weight  will  be  consumed  as  burned  fuel  in 
traveling  650  nautical  miles.)  Thus  the  specific  power  and  the  load- 
carrying capacity  of  the  vehicle  completely  determine  the  range-payload 
characteristics  and  direct  operating  costs. 

We  will  consider  the  problem  of  estimating  the  propulsion  specific 
power  in  later  Sections. 

NON-UNIFORM  NOZZLE  PARAMETERS 

In  the  derivation  of  the  cushion  performance  relationships 


(Equations  [B-4  ] through  [B-10]),  it  was  implied  that  the  quantities 
t,  8,  h,  and  p were  uniform  around  the  periphery. 

NOZZLE  THICKNESS  PARAMETER  x UNIFORM  - Actually  the  results  are 
valid  under  a somewhat  less  limiting  restriction;  namely,  the  restric- 
tion that  the  parameter 

x = (1  + sin  0) 

h 

be  uniform  around  all  parts  of  the  free  periphery,  and  that  air  be 
supplied  to  all  parts  of  the  nozzle  at  uniform  total  pressure  pt- 
However,  if  the  jet  angle  0 varies  around  the  periphery,  then 
it  is  necessary  to  use  an  average  value  of  the  discharge  coefficient, 
calculated  as  follows:  We  define  a curvilinear  coordinate  c,  measured 

along  the  cushion  periphery  (in  plan  view),  such  that 

£ dc  S'  C 

j h dc  = Sg 

Then  we  define 

* 

(P)  si-  £ P h dc 
v c'av  S J c 
g 

and  insert  this  value  in  place  of  in  Equations  [b-5]  through  [B-10]. 

(In  other  words,  we  take  a weighted  average  of  £ around  the  periphery, 
with  air  gap  h as  weighting  factor.) 

PRODUCT  £p  ( 1 - e ^X)j  UNIFORM  - A still  further  generalization  is 
possible.  If  the  air  is  supplied  to  different  parts  of  the  nozzle  at 
different  total  pressures  p (that  is,  from  different  compressor  systems 
all  operating  at  the  same  compressor  efficiency  and  duct  efficiency  Ty) , 

Equations  [B— 4 ] through  [B-10]  are  still  applicable  provided  the  quantity 
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is  uniform  around  all  parts  of  the  free  periphery.  Furthermore,  it 
is  assumed  that  the  above  weighted  average  will  be  used  in  place  of 
pc>  and  a weighted  average 


\ 


J5 

Ap/p 


/ — c 


fill  be  used  in  place  of 


t 

0C 
Ap/ p' 


av 


■ f f 


A?7iT 


h dc 


This  generalization  is  of  practical  importance,  for  example,  in 
case  peripheral  jets  are  to  be  used  around  most  of  the  free  periphery, 
while  the  cushion  is  allowed  to  function  as  a plenum  type  at  the  stern. 
(This  is  one  way  of  avoiding  the  tendency  of  stern  nozzles  to  snag  on 
obstacles  or  scoop  up  water.)  Equations  [B-4j  through  [B-10]  are  appli- 
cable to  this  case,  using  the  weighted-average  parameters  ^ and 

. c j av 

( jyc  \ 

l ~Ap7p  ) ’ Prov^ed  the  ait  which  escapes  under  the  stern  is  con- 

tinually  replaced  by  a separate  compressor  delivering  air  directly  to 
the  cushion  at  a total  pressure  pf  = Ap. 


OFF-DESIGN  CUSHION  PERFORMANCE 

When  the  conditions  specified  in  the  preceding  paragraphs  are  met; 

l -2x 

that  is,  either  x and  p^  are  uniform  or  Pt\l  “ e j is  uniform  around 
the  free  periphery,  the  cushion  is  said  to  be  operating  "on-dcsign." 
Roughly  speaking,  the  physical  significance  of  on-design  operation  is 
that  the  nozzle  conditions  at  each  section  of  the  periphery  are  just 
sufficient  to  maintain  the  cushion  pressure  Ap  without  having  air  entering 
or  leaving  the  cushion  itself. 

When  these  conditions  are  not  met  (off-design  cushion  operation), 
the  jet  is  relatively  stronger  in  some  sections  and  weaker  in  others; 
and  equilibrium  is  established  through  the  mechanism  of  air  entering 
the  cushion  from  the  strong  sections  of  the  jet  and  escaping  under  the 
weak  sections.  (Note  that  wo  define  off-design  cushion  operation  inde- 
pendently of  whether  or  not  the  compressors  are  operating  on-design. 
Compressors  will  be  discussed  in  the  next  section.) 
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Certain  conclusions  can  be  drawn  immediately. 

1.  Considering  peripheral-jet  machines  in  operation  over  a level 
surface : 


a.  If  h,  t,  and  9 are  uniform  around  the  free  periphery,  then 
the  cushion  will  operate  on-design  so  long  as  it  I3  in  level  attitude, 
regardless  of  operating  height. 

b.  If  h,  t,  or  9 varies  around  the  free  periphery,  then  only 
one  operating  height  and  attitude  can  provide  on-design  cushion  operation. 
(There  are  mathematical  exceptions  to  this  rule,  hut  not  practical  ex- 
ceptions.) 


c.  Practically  speaking,  per ipheral- jet  machines  cannot  provide 
on-design  operation  over  an  irregular  surface. 


2.  In  principle  (and  In  practice,  within  reason),  plenum  type 
machines  provide  on-des ign  cushion  operat  ion  regardless  of  operating 
height , attitude,  or  surface  irregularities . 

In  other  words,  peripheral-jet  cushions  practically  never  operate 
on-design,  whereas  plenum  cushions  practically  always  operate  on-design. 
We  must  therefore  re-examine  our  previous  conclusion,  that  peripheral- 
jet  types  are  more  efficient,  in  light  of  the  penalty  to  be  paid  in 
practice  due  to  off-design  conditions. 

It  is  fairly  easy  to  see  that  there  is  a penalty  for  off-design 
operation.  It  is  rather  difficult  to  calculate  the  precise  magnitude 
of  the  penalty.  We  will  consider  only  a few  of  the  simpler  cases,  In 
order  to  get  some  idea  of  the  magnitude  of  the  penalty. 

TWO-DIMENSIONAL  GEM  IN  PITCH  - When  a sidewall  craft,  on-design  J « 
level  trim,  is  placed  at  a finite  trim  angle,  a part  of  the  air  from  the 
low-end  nozzle  will  flow  through  the  cushion  and  pass  out  under  the  high- 
end  jet. 


If  the  cushion  pressure  Ap,  total  pressure  pt>  and  volume  flow 
rate  Q remain  unchanged,  then  the  relationship  between  the  air-gap 

decrease  Ah^~^  at  the  low  end  and  Increase  Ah^  at  the  high  end  can  be 
estimated,  In  the  case  x « 1,  as  follows: 

The  momentum  balance  across  the  low-end  .jet  gives 


Ap  ( h - Ah ^ = Ap  h - 2 p Q1  /-  p 

' v p c 


4„<->  iid'/t  p 

Ap  V P t 


x « 1 


If  we  assume  the  total  pressure  of  the  portion  Q'  of  the  air  is  reduced 
to  Ap  in  traversing  the  cushion,  we  have: 

J\  * 


(-)  S 

to—'  = 4 /il 

Ah(t)  J to 


or,  since  Ap  = 2x  pfc  when  x « 1, 

a£\  J- 


This  relationship  Is  plotted,  along  with  a modified  relationship  giving 

h<“> 

the  correct  behavior  — — 1.0,  as  x -• 

h''  ; 


1 1 r — t — r ~t~i  n 

Two-Dlraens lonal  GEM 
■v  Height  Loss  Due  to 
Pitch 


.tanh  ^ 


.4  .6  .8  j 


x - £ (1  + tin  8) 


2 3 4 


We  cannot  rely  very  heavily  on  either  of  the  relationships  pre- 
sented in  the  graph,  except  in  cases  of  rather  thin  nozzles  (small  x); 
but  evidently  the  loss  of  average  air  gap  with  finite  pitch  angles  is 
substantial  in  such  cases.  (It  is  interesting  to  note  that  this  behavior 
would  give  rise  to  a strong  coupling  between  pitch  and  heave,  in  case 
of  dynamic  motions;  that  is,  a pitching  oscillation  would  excite  heave 
oscillations.) 

OBSTACLE  CLIMBING  - Consider  a sidewall  craft  attempting  to  tra- 
verse an  obstacle  just  the  proper  width  to  pas3  between  the  sidewalls. 


For  simplicity  we  assume  that  level  trim  is  maintained  and  total  pres- 
sure p remains  constant  as  the  obstacle  is  approached.  In  the  situation 
sketched,  we  liave 


> x « 1 


= * (1  + sin  e) 
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I 


I 

I 


1 


: 


(Here,  h is  air  gap  corresponding  to  level-trim  operation  at  a distance 
from  the  obstacle.)  This  relationship  is  shown  graphically  below,  to- 
gether with  an  arbitrarily  modified  relationship  consistent  with  the 
obvious  limiting  case,  ^ - 1,  as  x -*  * 


.1  .2  .U  .6  .8  1 2 3 4 

x * ^ (1  + sin  9) 

Two-Dimensional  GEM  Step  Obstacle  Negotiable  In  Level  Trim 


FLIGHT  OVER  WAVES  - A problem  of  considerably  greater  practical  in-~ 

terest  is  posed  by  the  case  of  a full-peripheral- jet  craft  (no  sidewalls, 

etc.)  moving  without  water  contact  normal  to  the  wave  crests  of  an  ideal- 

Z Z 

ized  sinusoidal  sea.  We  consider  the  case  t-  » l amd  — » 1,  so  that  the 

b A. 

area  of  the  air  gap  is 


Sg  * 2 l (h  + 4h  - | ) 


Unfortunately,  as  is  often  the  case  when  we  come  to  a problem  of  prac- 
tical importance,  this  problem  is  considerably  more  complicated  than 
the  artificial  examples  of  the  preceding  two  paragraphs.  We  will  content 
ourselves  with  some  deductions  concerning  the  qualitative  behavior.  We 
will,  as  before,  adopt  the  reasonable  assumption  that  the  quantities 
Q and  pt  remain  fixed  (that  is,  independent  of  wave  height)  at  values 
producing  an  air  gap  of  h feet  in  level  trim  over  a level  surface. 

For  the  plenum,  x -•  ”,  it  is  clear  that 


Iherc  fore 


but 


S = Constant 
g 


h + Ah  = — + Constant 


Ah  = 0 when  H = 0 


so 


Also,  for  a very  thin  jet,  x -*  0,  it  is  clear  that  Ah  = 0,  since  the 
portion  of  the  air  gap  corresponding  to  Ap  in  the  wave  troughs  must  be 
filled  by  a volume  flow  of  air  diverted  from  the  jet  near  the  wave  crests 
and,  as  x -•  0 , the  volume  flow  of  air  available  also  approaches  zero. 
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This  behavior  la  graphed  below,  together  with  Intuitive  estimates  of  the 
behavior  when  x Is  finite.  It  must  be  borne  In  mind,  however,  that 
these  Intuitive  estimates  have  no  claim  to  validity,  quantitatively. 


PLENUM  VERSUS  PERIPHERAL  JET  - Our  considerations  of  off-design 
cushion  performance  have  raised  some  extremely  Important  questions  re- 
garding the  widely  accepted  notion  that  peripheral- jet  craft  are  more 
efficient  than  plenum  craft. 

There  Is  no  doubt  that  If  they  are  compared  on  the  basis  of  equal 
weight,  cushion  area,  and  daylight  gap  area,  and  also  on  the  basis  of 
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operation  over  a level  surface,  then  the  peripheral  Jet  is  substantially 
more  efficient.  Equations  [b-10]  and  [B-10*  ] show  the  cushion  power  to 
be  roughly  proportional  to  the  quantity  4-  (Ap/pt);  and  Figure  B-3 

shows  that  this  quantity  is  likely  to  be  only  70  to  80  percent  as  large, 
for  well  designed  peripheral  jets,  as  for  plenums.  (Moreover,  we  will 
see  in  the  following  sections  that  the  peripheral  jet  lends  itself  to 
slightly  smaller,  lighter  compressors  than  the  plenum,  and  produces  less 
momentum  drag.) 

However,  comparisons  on  the  bas is  of  operation  over  level  surface 
with  equal  daylight  gap  area  are  fallacious  1 The  reason,  of  course,  is 

that  a significant  daylight  gap  is  of  no  value  in  operation  over  a level 
surface.  The  value  of  the  daylight  gap  1 ies  entirely  in  avoiding  or 
minixlng  the  contact  with  an  irregular  surface. 

A more  proper  comparison  should  be  made  on  the  basis  of  operation, 
with  equal  daylight  gap  area,  over  a representative  irregular  surface. 
Unfortunately  this  is  not  easy  to  do.  Even  when  we  chose  very  simple 
irregular  surfaces  for  the  examples  of  the  preceding  paragraphs,  we 
did  not  succeed  in  determining  the  off-design  cushion  performance  for 
peripheral  jets  with  sufficient  accuracy  to  support  anything  more  than 
very  general  conclusions.  What  we  did  learn  can  be  summarized  as  follows: 
At  constant  power,  the  plenum  (x  °°)  provided  the  same  daylight  gap 
area  over  an  irregular  surface  as  over  a level  surface;  whereas  the  thin 
peripheral  jet  (x  — 0)  lost  as  much  as  half  its  daylight  gap  area  in 
moving  from  a level  to  an  irregular  surface.  (In  fact,  we  could  easily 
construct  other  examples  in  which  the  loss  would  be  more  than  half.) 
However,  the  crucial  question,  of  how  much  is  lost  with  intermediate 
thicknesses,  of  the  order  of  x = 1,  say,  remains  inadequately  answered. 
Almost  certainly,  a noticeable  fraction  — and  perhaps  most  — of  the 
peripheral  jet's  apparent  20  to  30  percent  level-surface  advantage  is 
bound  to  disappear  when  one  succeeds  in  making  proper  comparisons;  and  the 
'’optimum"  nozzle  thickness  parameter  will  undoubtedly  be  found  to  be 
still  thicker  than  the  value  of  x = 0.7,  which  would  be  deduced  from  Fig- 
ure B-3. 


Thus,  the  question  of  which  is  better,  plenum  or  peripheral  Jet, 
may  be  largely  an  academic  question.  If  the  answer  is  peripheral  Jet, 
it  will  likely  be  a very  plenum-like  peripheral  jet,  with  performance 
very  like  a simple  plenum. 

It  is  extremely  important  to  bear  in  ra^nd,  however,  that  peripheral 
Jets  with  thin  nozzles  will  suffer  a severe  penalty  in  off-design  cushion 
operation.  The  importance  of  keeping  this  in  mind  is  emphasized  by  the 
facts  that  we  have  adopted  cushion  performance  formulas,  Equations  [B-4  ] 
through  [b-10],  which  are  based  upon  on-design  operation;  r^d  that  con- 
siderations of  the  internal  aerodynamics  and  the  momentum  drag  provide 
arguments  in  favor  of  thinner  nozzles,  as  we  shall  see  in  the  following 
sections . 

Finally,  it  is  necessary  to  remember  that  the  whole  concept  of  on- 
design  cushion  operation  (and  on-design  cushion  performance  as  given  by 
Equations  [B-4 ] through  [B-10])  is  only  a convenient  idealization.  In 
actual  practice,  air  cushions  (except  for  simple  plenum  types)  practically 
never  operate  on-design. 
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Cushion  Pressure  Ratio 


Figure  B-l  - Air  Cushion  Pressure  Ratio  as  a Function  of  Nozzle 

Thickness  Parameter 
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Parameter  and  Peripheral  Jet  Efflux  Angle 
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Foreword 


This  report  is  based  on  a lecture  series  presented  by  the 
I authors  at  the  von  Karman  Institute  for  Fluid  Dynamics,  Rhode- 

Saint-Genese , Belgium,  in  May  1965;  and  at  the  University  of 

i 

Maryland,  College  Park,  Maryland,  in  July  1965.  The  lectures  were 


prepared  under  the  joint  auspices  of  the  David  Taylor  Model  Basin 
and  the  Naval  Air  Development  Center.  They  were  presented  in 
Belgium  under  the  joint  sponsorship  of  the  von  Karman  Institute 
and  the  Advisory  Group  for  Aerospace  Research  and  Development 
(AGARD) ; and  in  Maryland  under  the  sponsorship  of  the  Assistant 
Secretary  of  the  Navy  for  Research  and  Development. 

The  revised  lectures  will  be  presented  as  follows: 

A.  Introductory  Survey 

B.  Air  Cushion  Mechanics 

C.  Internal  Aerodynamics 


F.  Cushion  Contributions  to  Stability 

G.  Seakeeping 

H.  Performance  Summary 


NOTATION 


R compressor  reference  dimension:  effective  radius  from 

the  compressor  axis  to  the  compressor  blades,  ft 

Q compressor  rotational  speed,  radians/sec 

V,  reference  blade  velocity  (OR),  ft/sec 

□ 

q^  reference  blade  dynamic  pressure  ( pa  V^3/2),  lb/ ft2 

V mean  radial  velocity  through  blade  row,  ft/sec 

n 

\ internal  advance  ratio  (V  /V.  ) 

n b 

compressor  flow  area  at  inlet  face,  ttJ 

Ss  compressor  reference  flow  area  at  blade  row  (including 

area  occupied  by  blades),  ft2 

S3  compressor  exit  area,  ft3 

S clishion  area,  ft3,  measured,  in  plan  view,  to  outer  edge 

of  nozzle  exit 


S daylight  gap  area  (cushion  perimeter  times  mean  daylight 

8 gap  (hC),  ft3) 

p air  density,  slugs/ft3 

a 

H.  . mean  inlet  total  pressure  loss,  lb/ft3 

*» 1 


reference  inlet  total  pressure  loss,  lb/ft2. 

Value  of  H.  . when  V =0 
Z,  i o 

duct  loss,  lb/ft2.  Mean  loss  of  total  pressure  between 
compressor  exit  and  nozzle  exit  (or  cushion  entry) 


'll 


qs 


compressor  inlet  reference  dynamic  pressure,  lb/ft2 


compressor  exit  reference  dynamic  pressure,  lb/ft2 


cushion  reference  dynamic  pressure,  lb/ft2 


m 


I,  t 


'l,  d 


inlet  loss  coefficient 


duct  loss  coefficient  ( ^/q3 


\ 


rr,  { /qo  H^>i 

ram  recovery  efficiency  TT 


N 

Z 


duct  efficiency 


\Apt  + \%J 

number  of  identical  compressors  in  operation 
number  of  compressor  blades 
blade  chord,  ft 


c 

AP 

Apt 

Q 

\ 

T1 


int 


jy 


compressor  solidity  ratio 


is 

2n  R 


blade  angle,  deg 

reference  relative  flow  angle  tan-1 
airfoil  angle  of  attack*  deg 
airfoil  sectional  lift  coefficient 

airfoil  sectional  drag  coefficient 

shaft  power  supplied  to  compressor,  lb-ft/sec 
cushion  pressure,  lb/ft3 

total  pressure  rise  across  compressor,  lb/ft3 
volume  flow  rate,  ft'1 /sec 


compressor  efficiency 


internal  efficiency 


/ Q Ap  \ 
l Pc  ) 

{\  \) 


cushion  discharge  coefficient 


iv 
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SUMMARY 

The  power  dissipated  through  losses  in  the  inlet-compressor-ducting 
system  is  discussed.  A simple  illustrative  derivation  of  the  character- 
istics of  a hypothetical  compressor  is  carried  out.  These  characteristics 
are  then  used  to  illustrate  the  construction  of  combined  internal  char- 
acteristics charts. 

Tiie  dominant  influence  of  the  air  gap  area  in  determining  the  neces- 
sary size  of  the  compressor (s)  is  emphasized.  . » 

INTRODUCTION 

The  "ideal"  cushion  power  is  the  product  of  the  volume  flow  rate 
times  the  total  pressure  at  nozzle  exit  (or  cushion  entry),  Q p . The 

actual  shaft  power  can  easily  be  several  times  larger  than  this  if  the 
compressor  and  ducting  are  not  properly  matched  to  the  cushion  require- 
ments . 

A detailed  cons ideration  of  compressor  design  and  performance  would 
be  beyond  the  scope  of  the  present  study.  On  the  other  hand,  some  facets 
of  the  internal  aerodynamics  of  GEM's  can  be  explored  most  conveniently 
in  terms  of  explicit  compressor  characteristics,  and  it  is  rather  important 
to  understand,  qualitatively  at  least,  why  the  compressor  behaves  as  it 
does.  To  this  end,  a simple  approximate  derivation  of  the  characteristics 
of  a typical  compressor  will  be  carried  out. 

ILLUSTRATIVE  COMPRESSOR  CHARACTERISTICS 

For  this  purpose  we  will  consider  only  one  type  of  compressor,  the 
airfoil-type  radial  compressor  sketched  in  Figure  C-l.  Moreover  we  will 
restrict  ourselves  to  the  case  of  moderate  solidity  and  moderate  internal 
advance  ratio  (Say,  c £ 0.6  X £ 0.4)  so  that  we  can  neglect  swirl,  mutual 
interference  between  airfoils,  and  higher  orders  of  X without  errors  so 
r,  ,;e  as  to  obscure  the  primary  effects. 

Within  these  limitations,  the  aerodynamic  radial  force  coefficient 
n blade  is  roughly  equal  to  the  sectional  lift  coefficient  of  the 
it  two-dimensional  airfoil  at  the  same  angle  of  attack  (see 
/ *« r*  C-l): 

= Radial  Force  per  Blade 
*"r  q,  (Blade  Area)  C l 


The  total  pressure  rise  across  the  blade  row  is 


(Radial  Force  per  Blade)  X (No.  of  Blades) 


APt  - Gcx  qb 


[C-l] 


Similarly,  the  tangential  force  coefficient  (see  vector  diagram, 
Figure  C-l)  is  given  approxunately  by: 


Tangential  Force  per  Blade 
(Blade  Area) 


X CI  + cd 


The  shaft  power  required  to  overcome  the  tangential  blade  forces  is: 
(Tangential  Force  per  Blade)  x (No.  of  Blades)  x V 

b 

- °<X  cl  + cd}  qb  Vb  Ss 

In  addition,  we  must  add  the  power  dissipated  in  friction  on  the 
housing  of  the  rotating  fan.  Taking  an  average  coefficient  of  friction 
we  can  express  the  power  dissipated  on  one  side  of  the  lower  disc 


1.26  R 


I pa  ( Vb  r)  Cf  2nr  dr  = 4Cf  q 


V R 
b b 


= 3Cf  qb  Vb  S* 

Using  arbitrary  estimates  for  C£  and  for  the  number  of  equivalent 
half-discs  to  represent  the  rotating  housing,  we  obtain  for  the  total 
shaft  power  input 


a c . q.  V,  S X + “ 
1 ’b  b s c 


[C-2  ] 


The  useful  power  output  Is  Apt  Q , where  Q , the  volume  rate  of  flow,  Is: 


Therefore 


* = Vn  S2  = Vb  S2  X 


Q c a CX  qb  Vb  S2  X 


The  compressor  efficiency  Is 


APt  Q 


Cl  ° Cl 


Collecting  the  various  results,  we  have 


X + -^  + ^ 

cx  acx 


[C-3] 


c . d . 0.02 

— = O c X + + 

Vb  qb  s2  l ac£ 


[C-2  ' ] 


*Pt  . 

- a c , 

qb  1 


[C-l'l 


where  Equations  [C-2  ,]  and  [c-l/]  have  been  simply  rewritten  from 
Equations  [C-2]  and  [C-l]  in  a form  convenient  for  construction  of 
a nondimens ional  compressor  characteristics  chart.  The  construction 


i 


f 


I* 

» 

* 

f 


of  the  chart  proceeds 
value  of  o is  chosen; 
culated  from 


as  follows: 


c . and  c ,/c  . 
I d l 


For 

are 


selected  values  of  0 and 
read  from  Figure  02;  X 


o , a 
is  cal- 


X = tan"1  (0  - a) 

and  Equations  [C-l,])  [02'J,  and  LC-3  ] are  applied  to  obtain  one  point 
on  each  of  the  graphs  of  Figure  03.  This  process  is  repeated  for 
successively  different  choices  of  a until  complete  lines  are  traced  out. 

The  following  points  must  be  emphas ized  in  regard  to  the  approxi- 
mate theory  outlined  above,  and  the  compressor  characteristics  chart 
(Figure  03)  derived  therefrom: 


1.  The  compressor  which  we  have  considered  (Figure  C-l)  is  not  a 
real  compressor,  in  the  sense  of  having  ever  been  constructed  and  tested. 

It  is  a purely  hypothetical  compressor  chosen  for  an  illustrative  example 
(a)  because  it  lends  itself  to  a very  simple  approximate  theoretical 
analysis  and  (b)  because  the  theoretical  characteristics  so  derived  are 
quite  typical  of  the  characteristics  of  the  real  compressors  likely  to  be 
used  in  GEM  app 1 icat ions . (The  main  qualitative  differences  which  one 
would  probably  find  if  he  constructed  and  tested  our  hypothetical  compres- 
sor are  these:  The  pressure-rise  chart  of  Figure  C-3  would  show  a slightly 

higher  peak  in  the  curve  of  constant  0 for  each  successively  higher  value 
of  0,  due  to  effects  of  higher  orders  of  X,  which  we  have  negLected;  and 
the  power-input  chart  would  show  somewhat  lower  ordinates  in  the  region 

X ••’0,  due  to  swirl  effects,  which  we  have  neg<®bcted.) 

» 


2.  The  theoretical  method  of  analysis  we  have  used  is  not  suitable 
for  serious  design  application.  Far  more  accurate  (but,  necessarily, 
far  more  complex)  methods  are  readily  available  in  innumerable  textbooks, 
reports,  etc.  The  value  of  the  derivation  presented  here  is  only  that 
it  leads  to  predictions  of  the  compressor  characteristics  which  are  quali- 
tatively correct  and  which  are  comparatively  very  easy  to  understand. 
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DUCT  LOSSES 

Now  that  we  have  settled  on  a usable  description  of  the  compressor 
characteristics,  the  next  step  is  to  examine  how  these  characteristics 

■ 

must  be  matched  to  the  internal  flow  characteristics  of  the  vehicle  as  a 
whole . 

Consider  a schematic  diagram  of  the  whole  flow  system: 


V 

o 


•*>>-*■* 


Tabic  C-l 

Total  Gage  Pressure  and  Volume  Flow  Rate  at  Several  Stations 


Station 

Descript  ion 

' - 

Total  Cage  Pressure 

Volume  Flow  Rate 

(0) 

Free  Stream 

% 

- 

(1) 

Compressor 
Inlet  Face 

v H/,ia  \ v *Vi 

o 

Q 

(2) 

Blade  Row 

- 

Q 

(3) 

Compressor 

Exit 

q - H . .+  Ap 

'o  i9  L t 

Q 

(4) 

Nozzle 

q - H.  .+  Ap  - H.  p_ 

t £,d  rt 

Q 

(. 

"Daylight" 

Gap 

Q = oO  /-  Ap  S 
c •/  P 8 

(C) 

Cushion 

AP  - P^) 

— 

As  indicated  in  the  table,  it  is  convenient  to  separate  the  inlet 
loss  into  two  parts:  the  loss  which  occurs  with  zero  free-stream  velocitv, 

and  an  additional  loss  associated  with  free-stream  velocity: 


where  Tl  is  the  "ram  recovery  efficiency." 


Next,  it  is  convenient  to  define  loss  coefficients,  in  terms  of 
reference  dynamic  pressures: 


- 1.7  R4*  N 


' [*  + rt  ^ © («i.„ + *•*  *1>t| 

Further,  assume  K = 0.02  and  K,  . ■ 1.0.  (The9e  assumptions  are 

i/|  1 ™ “ 

reasonable  for  a typical  peripheral-jet  design  with  ducting  which  con- 
sists of  nothing  more  elaborate  than  a large  unobstructed  "plenum" 
through  which  the  air  flows  to  reach  the  noc?.les;  and  also,  for  .a 
plenum-type  GEM  in  which  air  is  delivered  more  or  less  directly  from  the 
compressor  to  the  cushion.  The  assumption  that  K - 1 is  approximately 

ij  > u 

equivalent  to  assuming  that  all  of  the  kinetic  energy  in  the  compressor 
efflux  is  lost  through  viscous  action.)  Then, 


1 + ^2.  a 


It  is  possible  to  understand,  from  this  formula,  one  of  the  funda- 
mental facts  about  GEM's  (which,  by  the  way,  has  led  to  the  downfall  of 
many  an  amateur  builder):  In  moderate-sized  GEM's  with  significant  day- 

light clearance,  the  compressors  require  a very  great  deal  of  space . 

Let  us  add  some  more  specifics  to  our  example:  Suppose  = 0.8, 

Pt 

& = 0.5,  the  average  daylight  clearance  (in  feet)  is  h,  the  cushion  is 
rectangular  with  length/bcam  ratio  I/b  = 2,  and  there  is  only  a single 
compressor  (N  = 1). 


Then  we  can  write 


^ = 


1 + 2.74 


1.235  /bh  ( ~ 


Tabulating  from  this  formula,  approximately,  for  h = 0.5  foot, 
the  maximum  diameter  cf  the  compressor  (2.52  R)  is  given  for  various 
combinations  of  beam  and  duct  efficiency  . 


It  wi.ll  be  very  much  worth  while  to  spend  a few  minutes  contem- 
plating this  table  and  the  development  which  led  up  to  it.  Of  the 
various  assumptions  which  were  U3ed  in  constructing  Table  C-2,  only  two 
can  be  manipulated  over  a sufficiently  wide  range,  in  a practical  design, 
to  really  strongly  influence  the  result.  These  are  the  assumptions  that 

h = 0.5  foot  and  K.  , = 1.0. 

2,d 

The  choice  of  daylight  clearance  h is  governed  primarily  by  con- 
siderations of  tradeoff  between  cushion  power  and  propulsion  power  (to 
be  discussed  in  the  succeeding  sections);  but  obviously  one  cannot  afford 
to  ignore  the  consideration  of  space  occupied  by  the  compressor.  Observe 
that  the  compressor  dimensions  vary  in  proportion  to  /h"  . That  is,  if  h 
is  increased  from  0.5  foot  to  2.0  feet,  the  compressor  diameters  in 
Table  C-2  mist  be  doubled;  and  if  h is  decreased  from  0.5  foot  to  0.125 
foot,  the  diameters  can  be  halved. 

Similarly,  the  compressor  dimensions  vary  approximately  in  propor- 

tion  to  ( K + 0.09)  . For  example,  if  the  duct  loss  coefficient 

\ / 

K a , could  be  reduced  from  1.0  to  0.23,  the  compressor  diameters  could 

*,d 

be  reduced  to  three-fourths  of  the  values  listed  in  Table  C-2.  This  is 
not,  by  any  means,  impossible;  but  it  nearly  always  turns  out  to  be 
impractical  in  actual  designs.  It  requires  smooth  ducts,  with  no  sharp 
bends  or  sudden  changes  in  stream-tube  area,  leading  from  the  compressors 
to  the  nozzles.  The  difficulty  is  that  to  do  this  usually  requires 
placing  the  compressors  in  locations  which  interfere  with  efficient  uti- 
lization of  space  for  payload.  For  example,  for  a single-compressor 
design,  highly  efficient  ducting  is  very  difficult  unless  the  compressor 
is  placed  near  the  center  of  the  planform.  It  is  sometimes  said  that 
the  need  for  efficient  ducting  can  be  avoided  by  the  use  of  an  efficient 
diffuser  to  decelerate  the  compressor  efflux,  but  this  is  deceptive. 

When  one  considers  the  total  space  occupied  by  the  compressor  and  diffuser 
combined,  It  usually  turns  out  to  be  more  economical  of  space  to  use  a 
larger  compressor  with  no  diffuser.  (The  same  conclusion  does  not  neces- 
sarily apply  to  considerations  of  weight  and  economy.) 
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Of  course,  each  design  presents  a slightly  different  problem,  and 
one  should  not  be  too  ha3ty  In  assuming  that  sophisticated  ducting  and/or 
diffusers  will  not  be  worth  while.  The  point  is  that  it  is  far  more 
dangerous  to  naively  assume  that  a value  of  K substantially  less  than 

unity  will  be  achieved,  unless  and  until  one  has  very  carefully  deter- 
mined precisely  how  this  will  be  achieved  and  what  it  will  cost.  In 

fact,  the  careless  designer  may  well  end  up  with  a value  of  K,  much 

z,  d 

greater  than  unity,  if  he  fails  to  assure  that  there  is  adequate  space 
in  his  ducting  for  the  stream  tubes  to  reach  their  destination  wit:. out 
"squeezing  through"  constricted  passages  before  reaching  the  nozzle  (in 
the  case  of  peripheral-jet  types)  or  the  air  gap  (in  the  case  of  plenum 
types) . 

Finally,  let  us  reflect  for  a moment  on  the  implications  of  Table 
C-2,  if,  instead  of  a single  large  compressor,  one  uses  a number  of 
smaller  compressors.  In  fact,  the  designer  usually  will  choose  a multiple- 
compressor  design,  except  in  cases  of  very  small  vehicles,  for  a variety 
of  reasons:  weight  and  balance,  reliability,  economy,  convenience  of 

fabrication  and  maintenance,  power  transmission  problems,  engine-compressor 
matching,  etc.  However,  use  of  multiple  compressors  does  not  generally 
alleviate  the  problem  of  space  occupied  by  the  compressors;  just  the 
contrary,  in  fact.  For  example,  if  four  compressors  were  used,  the  diam- 
eter of  each  would  be  half  the  value  listed  in  Table  C-2.  Taking  into 
account  that  these  four  half-size  compressors  must  be  spaced  far  enough 
apart  to  avoid  excessive  interference  between  their  respective  effluxes, 
it  is  easy  to  see  that  they  will  generally  entail  more  compromise  of  the 
useful  space  aboard  the  vehicle  than  one  full-size  compressor  would  have 
done . 


The  following  conclusions  are  Indicated: 


1.  The  duct  efficiency  depends  on  a number  of  design  details, 
but  is  by  far  most  strongly  dependent  on  the  ratio  of  "daylight"  gap 
area  to  compressor  exit  area,  Sg/S3  . 


2.  The  proportionate  planform  space  occupied  by  compressors  (if 

good  duct  efficiency  is  maintained)  is  thus  primarily  dependent  on  the 

ratio  S /S  . Roughly  speaking,  this  space  problem  tends  to  become 

oppressive  when  S /S  > 0.05,  say,  unless  unusual  measures  are  taken  to 
6 

avoid  duct  losses. 


OVER-ALL  INTERNAL  EFFICIENCY 

Equation  [c-4  ] expresses  the  duct  efficiency  in  terms  of  relation- 
ships between  the  dimensions  of  the  compressor  and  the  dimensions  of 
the  vehicle  as  a whole.  It  will  now  be  instructive  to  examine  the  over- 
all internal  efficiency. 

\nt  = \ \ 

We  have  an  approximate  expression  for  (Equation  [C-3])  for  a compressor 
of  the  type  shown  in  Figure  C-l.  It  will  be  convenient  to  re-express 
Equation  [C-4 ] in  terms  of  the  compressor  performance  parameters.  This 
can  be  approached  directly,  but  it  is  considerably  easier  to  start  over 
from  the  beginning.  ' 

Let  us  write 


but 


Apt-  V “ %,d+  \qo 

° > 

Apt  " km  qi  - + ^ q0 


q 


3 


q 

i 


x.  q,. 
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71 . = 1 - 0.689  - — 7~ 

d Apt/ql 


which  allows  us  to  assign  a value  of  T),  immediately  to  each  point  of 

d 


the  pressure-rise  chart  (Figure  C-3)  and  to  plot  the  product  Tj^  T|  = T)  , 
instead  of  , on  the  efficiency  chart. 

Next,  we  can  convert  the  pressure-rise  chart  of  Figure  C-3  to 
a cushion- pressure  chart,  and  thus  obtain  a system  characteristics 
chart  for  the  whole  vehicle,  analogous  to  the  compressor  character- 
istics chart  for  the  compressor  alone.  This  conversion  proceeds  as 
follows.  Note  that 
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where  the  subscript  "demand"  denotes  that  this  relationship  is  derived 

from  the  fact  that  a certain  cushion  pressure  Ap  acting  across  gap  area 

S with  discharge  coefficient  <£f  demands  a corresponding  value  of  X V.  S , 
g c b a 

from  continuity  considerations.  Wc  need  also  an  "availability"  relation- 
ship, giving  the  cushion  pressure  which  the  compressor  can  sustain.  This 
is  easily  obtained  as  follows 


Ap 


q 


b 


or 


Ap  \ Ap  Apt 

qb  / available  Pt  qb 


\ %/%  \ 
Apt/qb  J 


"Internal  System  Characteristics  Charts"  constructed  from  these 
relationships  are  presented  in  Figure  C-4.  These  charts  provide  fairly 
complete  information  on  the  internal  system  performance,  and  inter- 
relationships between  the  main  design  parameters,  subject  to  the  simpli- 
fying assumptions  which  have  been  introduced. 

The  following  points  are  worth  noting,  in  regard  to  Figure  C-4: 

1.  Within  the  range  considered,  most  aspects  of  internal  per- 
formance improve  with  increasing  solidity,  o.  We  have,  therefore,  not 
established  an  "optimum"  solidity.  We  may  surmise,  however,  that  the 
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performance  shown  for  o “ 0.6  is  not  significantly  different  from 
optimum.  Higher  solidities  were  not  investigated  because  a more  complex 
compressor  theory  would  have  been  required. 

2.  For  o = 0.6,  there  is  clearly  an  "optimum"  combination  of  blade 
angle  3 and  advance  ratio  X , from  the  standpoint  of  providing  maximum 
internal  efficiency.  Figure  C-4  does  not  include  enough  values  of  p 

to  define  this  optimum  combination  precisely,  but  it  is  near  ft  - 20°, 

X = 0.22. 

3.  One  usually  should  not  choose  this  optimum  combination  as  the 
design  point,  because: 

(a)  The  optimum  point  lies  uncomfortably  close  to  the  stall 

point,  and  intermittent  stall  might  be  experienced  in  unsteady  operating 
conditions.  (Roughly  speaking,  for  a given  p the  compressor  is  stalled 
at  all  values  of  X lower  than  the  peak  of  the  curve,  and  is  said  to 
be  "full  stalled"  when  X = 0.)  ^ 

(b)  If  the  design  point  corresponds  too  closely  to  the  maximum 
pressure  output  of  the  compressor,  the  full-stalled  pressure  output 
might  be  inadequate  for  take-off. 

(c)  This  choice  of  the  design  point  might  preclude  operating 
the  lift  system  at  less  than  full  power,  or  preclude  flight  at  gross 
weight  greater  than  the  design  value  (most  desirable  capabilities,  because 
they  provide  extended  range  and/or  payload  in  unusually  favorable  oper- 
ating conditions).  (Note,  in  this  regard,  the  tremendous  advantage 
which  a variable-g  compressor  would  enjoy.'  Review,  also,  Figure  C-.3 

from  this  viewpoint.) 

(d)  Careful  study  of  the  relationships  illustrated  in  the  Cushion 
Pressure  Chart  of  Figure  C-4  will  reveal  that  choice  of  hlgher-than- 
optimum  values  of  both  g and  X permits  use  of  smaller  compressors,  with 
savings  in  weight  and  space.  The  tradeoff  of  internal  efficiency 

versus  compressor  weight  and  size  is  an  extremely  important  one.  We 
will  not  attempt  to  formulate  this  tradeoff  quantitatively  herein,  because 
the  quantitative  benefit  of  a given  weight  or  size  reduction  is  very 
difficult  to  determine,  and  it  varies  greatly  from  one  application  to 
another.  This  tradeoff  must  be  considered  with  utmost  care,  however, 
in  actual  design  efforts. 


At  the  rtsk  of  becoming  monotonous,  a final  reminder  must  be  given 
of  the  assumptions  and  simplifications  which  led  up  to  Figure  C-4: 

1.  Only  one  type  of  compressor  (Figure  C-l)  was  considered,  and 
only  limited  ranges  of  its  allowed  variables  (g,  X,  and  o)  were  con- 
sidered. A somewhat  more  complete  picture  could  be  obtained  by  con- 
sidering a broader  range  of  variables,  but  this  would  require  introduc- 
tion of  mote  complex  theory.  In  defense  of  the  presentation  given,  it 
may  be  pointed  out  that  the  system  characteristics  derived  herein  are 
qualitatively  representative  of  most  practical  systems.  In  fact,  the 
simple  theory  for  axial  compressors  is  exactly  analogous  to  that  given 
for  the  airfoil-type  radial  flow  compressor,  the  main  first-order 
difference  being  a different  numerical  relationship  between  the  areas 

S , S , and  S ; and  the  centrifugal  compressor  is  simply  an  extremely 
is  3 

high-solidity  version  of  the  compressor  considered  here. 

2.  Only  one  of  the  infinite  number  of  possible  combinations  of 

duct  loss  coefficients  is  reflected  in  Figure  C-4,  namely  K.  . = 0.02, 

t,  l 

K = 1.0.  This  combination  is  quite  representative  of  current  vehi- 

jCj  J Q 

cles,  but  it  must  be  remembered  that  current  vehicles  represent  a rela- 
tively primitive  state-of-the-art.  (Hovercraft  were  first  seriously 
considered  around  1957!)  In  a previous  discussion,  the  merits  of  more 
sophisticated  ducting  were  somewhat  discounted,  fr<jm  the  standpoint  of 
alleviating  the  space  problem  imposed  by  large  compressors.  However, 
the  same  arguments  do  not  necessarily  apply  to  questions  of  compressor 
weight  and  economy.  The  relationships  given  in  this  section  permit  the 
reader,  with  fairly  nominal  effort,  to  construct  other  "system  charac- 
teristics charts,"  similar  to  Figure  C-4,  but  with  different  assumptions 
regarding  duct  loss  coefficients;  and,  if  he  wishes,  he  could  use  manu- 
facturers' compressor  characteristics  in  place  of  the  illustrative  char- 
acteristics, derived  herein.  An  exercise  of  this  kind  is  definitely 
recommended  for  readers  with  serious  design  intentions. 
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Foreword 

This  report  is  based  on  a lecture  series  presented  by  the 
authors  at  the  von  Karman  Institute  for  Fluid  Dynamics,  Rhode- 
Saint-Cenese,  Belgium,  in  May  1965;  and  at  the  University  of 
Maryland,  College  Park,  Maryland,  in  July  1965.  Hie  lectures  were 
prepared  under  the  joint  auspices  of  the  David  Taylor  Model  Basin 
and  the  Naval  Air  Development  Center,,  They  were  presented  in 
Belgium  under  the  joint  sponsorship  of  the  von  Karman  Institute 
and  the  Advisory  Group  for  Aerospace  Research  and  Development 
(AGARD) ; and  in  Maryland  under  the  sponsorship  of  the  Assistant 
Secretary  of  the  Navy  for  Research  and  Development. 

The  revised  lectures  will  be  presented  as  follows: 


A. 

Introductory  Survey 

B. 

Air  Cushion  Mechanics 

C. 

Internal  Aerodynamics 

D. 

Drag 

r-  • 

F, 

Cushion  Contributions 

to  Stability 

G. 

Seakeeping 

H. 

Performance  Summary 

ll 


NOTATION 


1 


j 

wave  amplitude,  ft.  (Half  of  the  trough  to  crest  height.) 
beam  of  pressure  region  or  bubble,  ft 

external  aerodynamic  drag  coefficient  (based  on  cushion 
area,  S) 

trunk  drag  coefficient  (trunk  drag/qaS) 

skin  friction  coefficient  . • • 

lift  coefficient  (lift/qaS) 

external  aerodynamic  drag,  lb 

ram  or  momentum  drag,  lb 

additional  sidewall  drag,  lb 

secondary  wavemaking  drag  (drag  due  to  retaining  air  bubble),  lb 
trunk  (or  skirt)  drag,  lb 

total  drag,  lb.  (The  appropriate  sum  of  the  component  drags.) 
wavemaking  drag,  lb 

nozzle  discharge  coefficient  of  ACV  base 
length  Froude  number  {y  I vgi)  (nond imens ional) 
area  Froude  number  ( V /.Jg  ,/ST)  (nondimens ional) 

/ °w/L  \ 

wave  drag  parameter  |f^  = — J (see  Figure  D-2) 

wave  drag  parameter  (see  Figure  D-3) 

geometry  factor  (G  = 1 for  full-peripheral  ACV; 

G < 1 for  CAB  ACV) 

ill 


average  wave  height,  trough  to  crest,  ft 

height  or  daylight  clearance,  ft 

additional  sidewall  depth,  ft  (see  Figure  D-5) 

maximum  wave  depression  under  craft,  f'.  (see  Figures  D-5  and  D-6) 
lift,  lb  (L  = W) 

length  of  pressure  region  (or  bubble),  ft 

short  length,  ft 

sidewall  wetted  length,  ft 

shortened  wetted  length  of  ACV  with  daylight  clearance,  ft 
(see  Figure  D-6) 

number  of  sidewall  wetted  sides 
propulsive  power,  lb-ft/sec 

pressure,  lb/fta 

dynamic  air  pressure  Oi  P V2),  lb/ft2 

dynamic  water  pressure  (k  p V2),  lb/ft3 

w 

resultant  lift  force,  lb 

♦ 

cushion  base  area,  ft3 

area  of  discharge  (daylight  clearance  area),  ft2 
velocity,  ft/sec 

jet  velocity  of  air  discharge  from  ACV  base,  when  fully 
contracted  , ft/sec 

vehicle  weight,  lb 

specific  weight  of  ACV  ^ ^ » lb/ft2 


X 

pa 

y 


vehicle  pitch  angle  or  angle  of  attack,  deg 

propulsive  efficiency 

wave  length,  ft 

mass  density  of  air,  slugs/ft3 

mass  density  of  water,  slugs/ft3 

weight  density  of  water  (p^  g) , lb/ft3 


I 

I 
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SUMMARY 


Six  vehicle  drag  terms  are  derived.  These  apply  to  the  sidewall 
captured  air  bubble  air  cushion  vehicle  and  the  full-peripheral  air 
cushion  vehicle. 
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INTRODUCTION 

The  subject  of  the  drag  of  air  cushion  vehicles  is  broken  up 
into  six  parts.  The  various  components  of  drag,  which  apply  to  a 
skirted  full-pheripheral  air  cushion  vehicle  (ACV)  and  to  a captured 
air  bubble  (CAB)  type  of  sidewall  vehicle  with  only  slight  modifica- 
tion, are  as  follows:  wave  drag,  D ; a secondary  wavemaking  drag 

w 

associated  with  the  drag  of  that  amount  of  sidewall  required  to  retain 

the  bubble  (or  cushion),  D ; additional  sidewall  drag,  D ; the 

s jW  s » a 

trunk  (or  skirt)  drag  of  the  full-peripheral  ACV,  Dt  (a  geometry 

factor  of  less  than  one  is  used  for  the  fore  and  aft  ski  drag  in  the 

CAB  type  craft);  the  external  aerodynamic  drag  (air  profile  drag  made 

up  of  form  drag  and  skin  friction  drag),  D ; and  the  ACV  ram  or 

e 

momentum  drag,  . The  total  drag,  D^,  , is  an  appropriate  sum  of  the 
drags.  All  the  drag  terms  will  ultimately  be  referred  to  vehicle 
weight,  W (which  equals  the  vehicle  lift,  L) , to  give  the  drag/weight 
(or  drag/lift)  ratio. 

The  propulsive  power  is  the  product  of  the  total  drag  , the 
vehicle  velocity  V,  and  the  reciprocal  of  the  propulsive  efficiency, 

11^  . The  specific  propulsive  power  (nondimensional)  is  given  by: 
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WAVEMAKING  DRAG 

A region  of  pressure,  proceeding  at  some  speed  V over  a body  of 
otherwise  calm  water,  generates  waves.  If  the  pressurized  region  is 
considered  to  be  at  rest,  an'  the  otherwise  calm  water  has  a velocity 
of  -V,  then  the  wave  pattern  generated  is  fixed  in  space. 

Our  interest  in  this  wave  is  that  an  ACV  pressure  region  generates 

such  a wave,  and  further  that  the  pressure  region  ultimately  acts  back 

on  the  vehicle  to  give  a resultant  air  pressure  force  tilting  back 

somewhat  from  the  vertical.  This  tilted  force  has,  in  addition  to  a 

lift  component,  a drag  component,  which  will  be  called  a wavemaking 

drag,  D . 
w 

The  simplest  wavemaking  case  is  that  of  a two-dimensional  (infinite 
beam  and  infinite  aspect  ratio)  short  pressurized  region  of  length  A2, 
pressure  p,  and  velocity  V.  This  was  treated  by  Lamb  (Reference  1, 
pages  402-403)  and  by  Wagner  (Reference  2)  and  others.  Lamb  notes  that 
at  a distance  o-f  about  half  a wavelength  from  the  origin,  on  the  down- 
stream side,  a simple-harmonic  wave  profile  is  fully  established.  The 
amplitude  "Aa"  is  given  by 

Aa  = 2pA^  [D-i; 

P V2 

Note  that  the  amplitude  is  proportional  to  pA2>  and  to  the  reciprocal 

of  the  square  of  velocity.  The  wave  length  \ is  proportional  to  the 

square  of  velocity  as  follows: 

. 2tt  V2  r 

x-  — - »2- 

The  wavemaking  drag/lift  ratio  is  proportional  to  Aa/\,  and  is  given  by 

t 

D . 

= £_£_Ai 

L p V*  [D-3] 


Hence,  the  wavemaking  drag  is 


Pw  V* 


Note  that  the  wavemaking  drag  is  proportional  to  the  square  of  the 
pressure  p,  and  is  inversely  proportional  to  the  fourth  power  of  veloc- 
ity. On  the  logaritlimic  plots  (Figure  D-2,  for  example),  this  fourth- 
power  relationship  can  be  easily  verified  for  the  two-dimensional  case. 


1 


The  next  simplest  wavemaking  case  Is  that  of  a two  dimensional 
(infinite  beam  and  infinite  aspect  ratio)  uniform  pressurized  region  of 
finite  length  l . Lamb  (Reference  1,  p 404)  indicates  a method  for 
integrating  the  previous  results,  taking  into  account  not  only  ampli- 
tudes, but  also  phases  of  the  component  wqves  generated  by  different  and 
separated  A £ regions.  Lamb  says  that  "we  can  easily  deduce  the  requi- 
site formulae";  but  in  fact  the  process  is  rather  complex  and  tedious, 
and  even  the  final  formulas  are  long. 

It  does,  however,  yield  formulations  of  surface,  shap*  in  three 
regions,  ahead  of  the  pressure  region  of  length  £,  within  it,  and 
behind  it.  The  one  simply  stated  result  is  the  wavemaking  drag/lift 
ratio. 

r - i ■ i?i  (l  - cos  ld-4] 


where  h and  £ are  shown  in  Figure  D-l.  This  formula  reduces  to 

equation  [D-3]  when  £ is  very  small  (or  when  g 1 « v2»  of.  » l> 

8 — 


or 


>>  1).  Note  that,  since  lift  L is  proportional  to  pres- 


sure, the  wave  drag  D^  is  proportional  to  the  square  of  the  pressure 
p.  Figure  D-l,  taken  from  Lamb,  shows  the  surface  shape  taken  by 
the  water  for  £ equal,  to  — or  for  a Froude  number  V/^g  l of  unity. 
In  this  case  the  pressure  length  £ is  0.159  times  wave  length. 


The  most  difficult  of  the  wavemaking  problems,  that  due  to  a con- 
stant pressure  region  over  a finite  region,  is  the  most  useful.  Fortu- 
nately, this  work  has  been  accomplished  in  recent  years  at  the  David 
Taylor  Model  Basin  (Reference  3).  Figure  D-2  shows  these  results  for  a 
rectangular  distribution  of  constant  pressure.  The  wave  drag  parameter 
f^  plotted  on  the  ordinate  (drag/lift  ratio  divided  by  pressure/ length 
ratio  and  constants)  is  plotted  against  the  length  Froude  number  V/Jgt,  . 

It  can  be  noted  that  the  wave  drag  Dy  is  proportional  to  lift  L , to 

pressure  p,  the  wave  drag  parameter  f , and  is  reciprocally  proportional 

I 

to  pressure  (or  bubble)  length  £ . The  value  of  the  wave  drag  parameter 

B 

f depends  only  on  the  length/beam  (£/ b)  ratio  and  Froude  number  F . 

JL  JL 

The  speeds  of  maximum  f^  are  referred  to  as  hump  speeds  or  Froude  numbers. 
In  the  sub-hump  region,  the  slopes  of  the  lines  in  Figure  D-2  were 


I 
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adjusted  to  the  slopes  of  early  British  Hovercraft  experimental  data. 
(See  Reference  4.)  The  actual  below-hump  speed  shape  of  f^  Is  prob- 
ably different  from  that  shown  in  Figure  D-2.  The  above  simplification 
is  only  justified  by  the  lack  of  importance  of  the  below-hump  speed 
values  to  the  present  task  at  hand.  Secondary  drag  humps  and  valleys 
cai  exist , however , and  they  could  be  of  considerable  significance  for 
some  purposes. 

Figure  D-3  is  also  a plot  of  Newman's  results  (Reference  3)  for  a 
rectangular  pressure  region,  but  here  the  pressure  or  bubble  length  l 
is  replaced  by  ^'S  , where  S is  the  total  pressure  area.  This  replace- 
ment of  £ by  /S  is  done  both  in  the  ordinate,  abbreviated  fy§- , and 
the  abscissa  F , and  results  in  a numerically  different  plot  from 
that  of  Figure  D-2.  The  advantage  of  Figure  D-3  is  that  equal  area 
pressure  regions  S,  at  the  same  area  Froude  number  F /§' , have  equal 
speeds;  whereas , this  is  not  true  of  Figure  D-2.  This  could  prove  to 
be  of  some  advantage  in  the  comparison  of  vehicles  of  equal  planform 
areas  at  the  same  speeds. 

Figure  D-4  is  a plot  similar  to  that  of  Figure  D-3  but  for  ellip- 
tical planform  shapes. 

The  resultant  ratio  of  wave  drag  to  lift  (or  drag/weight)  is 
written,  based  on  Figure  D-2.  For  Figures  D-3  and  D-4,  similar  expres- 
sions apply.  , 


W 


-^7 

W1 
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If  external  aerodynamic  lift  C q S provides  partial  vehicle  support, 

Li  3 

then  the  cushion  pressure  is  less  than  the  vehicle  specific  weight  w, 
which  is  (W/S) . In  terms  of  the  vehicle  specific  weight  w,  the  wave 
drag/weight  ratio  is 


D 


w 
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SIDEWALL  DRAG  IN  A CAPTURED  AIR  BUBBLE  TYPE  OF  AIR  CUSHION  VEHICLE 

Figure  D-5  shows  a CAB  ACV  in  an  above-hump  speed  condition.  The 
angle  a is  exaggerated  to  point  out  two  regions  of  sidewall  wetting 
and  drag. 

A primary  function  of  sidewalls  (or  sideboards  or  side  skegs)  is 
to  prevent  loss  of  air  along  the  side  of  the  vehicle.  This  is  done 
by  a triangularly  wetted  sideboard  section  of  maximum  wetted  height  h^; 
that  is,  wetted  on  the  outside  only.  This  wetted  area  iS  referred  to 
as  sidewall  due  to  the  bubble.  The  assumption  here  s that  the  outside 
water  line  continues  horizontally;  whereas,  the  inside  water  line  is 
turned  down  by  the  angle  a. 

Below  this  triangular  region  of  sideboard  due  to  the  bubble,  is  a 
region  fully  wetted  on  both  sides  that  is  referred  to  as  additional 
sidewall , which  is  fully  wetted,  both  inside  and  outside. 

ADDITIONAL  SIDEWALL  DRAG  ~ The  amount  of  additional  sidewall 
depth  h (Figure  D-5)  in  smooth  water  would  be  held  to  a minimum  at 
high  speeds,  but  increased  for  transiting  the  hump  speed.  In  waves  at 
high  speed,  fore  and  aft  seals  or  skis  (skirts  or  trunks  at  zero 
"daylight"  clearance)  are  set  very  close  to  the  bottom  of  the  sidewalls, 
but  waves  of  height  H will  nevertheless  wet  an  equivalent  height  h 

a 

of  approximately  one-half  the  average  wave  height. 


h 
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The  reason  for  this  is  that  (for  the  above  ski  setting)  the  bottom  of 
the  sidewall  runs  at  the  wave  trough,  and  the  fore  and  aft  skis  or  seals 
move  up  and  down  rapidly,  staying  very  close  to  the  water  surface. 

The  additional  sidewall  drag  is  given  by  the  product  of  the  skin 
friction  coefficient  C^,  the  dynamic  pressure  in  water  q^,  and  the 
wetted  area.  Wetted  area  is  the  product  of  sidewall  length  Xg  , wetted 
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The  additional  sidewall 


height  h , and  number  of  wetted  sides  n. 
a 

drag/weight  ratio  Is  then,  for  n equal  4: 
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SECONDARY  WAVEMAKINC  DRAG  OR  SIDEWALL  DRAG  DUE  TO  RETAINING  THE 
AIR  BUBBIE  — Figure  D-5  shows  triangularly  wetted  sidewall  sections 
of  maximum  wetted  height  hfe  that  are  wetted  on  the  outside  only. 
These  triangular  sections  prevent  air  loss  along  the  vehicle  sides. 
Their  total  water  wetted  area  for  the  whole  vehicle  is  h^  • The 
drag  is  the  product  of  the  drag  coefficient  C^,  water  dynamic  pres- 
sure q , and  wetted  area  l h.  . 
w B b 
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From  Figure  D-5,  it  can  be  seen  that 


Thus  , 


W 


[D-10] 


D q.  D 

s ,w  B c w w 

W b r w W 


tD- 11] 


The  expression  (D-ll)  is  for  a CAB  type  of  ACV  with  sidewalls,  and 
no  appreciable  daylight  clearance  h. 
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Figure  D-6  shows  a sidewall  ACV  with  some  daylight  clearance  h, 
but  with  the  sidewall  still  wetted  to  some  extent  on  the  outside. 

The  wetted  height  is  (hj,  - h) ; the  wetted  length  is  l',  and 


[D-12] 


The  wetted  area  is  the  product  of  l'  and  (h^  - h} , 


This  equation  is  valid  for  a more  general  case  than  (D-ll).  It  is 
valid  for  a sidewall  ACV  with  some  daylight  clearance  h;  but,  never- 
theless, with  some  wetted  area,  as  shown  in  Figure  D-6.  When  h// 

B 

equals  and  then  exceeds  D^/W,  this  drag  becomes  and  remains  zero. 

SIDEWALL  WAVEMAKING  DRAG  - Because  of  the  high  fineness  ratios  of 
CAB  type  ACV  sidewalls,  and  low  immersed  depths  when  "on  the  bubble," 


sidewall  wavemaking  drag  Is  usually  negligible.  Nevertheless,  a form 
for  computing  It  Is  available,  based  on  wave-resistance  theory  for  thin 
bodies  (Reference  5).  The  resistance  R Is  given  in  terms  of  beam  B, 
draft  H,  length  L,  and  the  resistance  coefficient  r as  follows: 


R 


I 

TT 
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The  resistance  coefficient  r is  plotted  in  Figure  D-7  as  a function 
of  the  reciprocal  of  Froude  number  squared  for  a parabolic  thin  body 
form  and  a "full"  form. 

FORM  DRAG  — Form  drag  of  water-immersed  portions  of  sidewalls  is 
also  usually  negligible  because  of  the  high  fineness  ratios  and  low 
immersions.  The  profile  drag,  made  up  of  skin  friction  and  pressure 
form  drag,  then  approximates  the  skin  friction  drag.  One  form  of  drag 
estimate  by  Hughes  (Reference  6)  is: 


Pressure  Form  Drag  = 1-16I  Td-16] 

Profile  Drag  i 

where  t and  l are  body  thickness  and  length,  respectively. 

TRUNK  OR  SKIRT  DRAG  IN  FULL- PERIPHERAL  AIR  CUSHION  VEHICLES 

The  most  important  new  element  which  the  British  have  introduced 
in  air  cushion  vehicle  technology  is  flexible  trunks,  or  skirts,  as 
illustrated  in  Figure  D-8.  These  trunks  have  in  no  way  changed  the 
basic  principles  of  air  cushion  vehicles  but  they  have  had  a very  profound 
effect  on  tfie  design  choices  that  the  designer  is  likely  to  make.  The 
mest  obvious  effect  on  the  designer's  choice  is  in  the  question  of  the 
"daylight"  gap  (the  clearance  between  the  bottom  of  the  skirt  and  the 
water)  and  in  the  hard  structure  clearance  (clearance  between  the  hard 
structure  of  the  vehicle  and  the  water).  The  introduction  of  trunks 
has  caused  the  designer  to  choose  much  lower  values  of  the  daylight  gap 
and  much  higher  values  of  the  hard  structure  clearance.  These  choices 
involve  design  tradeoffs;  for  example,  the  choice  of  daylight  gap 


clearance  involves  the  tradeoff  between  lift  power  and  propulsion 
power.  Using  small  values  of  daylight  clearance,  the  lift  power  is 
greatly  reduced.  On  the  other  hand,  as  soon  as  the  vehicle  operates 
not  on  a smooth  surface,  but  in  a seaway,  the  flexible  trunk  drags 
through  the  water  and  there  is  an  increase  in  propulsive  power  re- 
quired . 

The  only  definitive  information  readily  available  on  the  subject 
of  the  drag  associated  with  pulling  flexible  trunks  through  waves  is 
represented  in  Figure  D-9.  It  comes  from  the  tests  of  the  VA-3 
Hovercraft  conducted  by  Republic  Aviation  under  contract  to  the  Office 
of  Naval  Research. 

In  the  left  graph  of  Figure  D-9,  a trunk  drag  coefficient,  re- 
ferred to  cushion  area  and  the  dynamic  pressure  of  air,  is  plotted 
versus  the  ratio  of  wave  height  to  cushion  length  for  two  different 
values  of  daylight  clearance.  It  was  then  surmised  that  air  cushion 
vehicles  tend  to  average  out  the  waves  over  which  they  fly,  and  there- 
fore significant  trunk  drag  should  not  occur  until  the  wave  height 
exceeds  twice  the  daylight  clearance.  On  this  basis,  in  the  right  graph 
in  Figure  D-9,  the  same  data  were  replotted  versus  the  ratio  of  the  wave 
height,  less  twice  the  daylight  clearance,  to  cushion  length.  It  was 
found  that  the  data  then  fell  nearly  on  a single  line,  which  could  be 
represented  by  a very  simple  nond imens ional  formula.  This  formula  is 
an  empirical  expression  of  the  VA-3  trunk  drag.  It  is  not  known  whether 
the  drag  characteristics  of  the  VA-3  trunks  are  favorable  or  unfavorable 
compared  with  other  designs.  Since,  however,  these  are  the  only  defin- 
itive data  available,  this  formula  will  be  applied  directly  to  estimate 
performance  in  a sea  state. 

The  trunk  drag  coefficient  , is  a function  of  the  average  wave 

height  H,  the  ACV  length  X,  and  the  daylight  clearance  h.  The  drag 
coefficient  is  the  drag  per  unit  cushion  area  S,  and  per  unit  air 

dynamic  pressure  q&. 
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The  trunk  drag/weight  ratio  is  then. 
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CAPTURED  AIR  BUBBLE  AIR  CUSHION  VEHICLE  OR  SKI  DRAG 

The  fore  and  aft  seals  or  skis  on  CAB  ACV's  and  models  have  been 
of  a mechanical  type  principally,  but  the  advantages  of  low-weight 
fabric  materials,  relative  to  high-frequency  response,  are  evident. 

For  this  reason  an  analogy  can  be  drawn  to  full-peripheral  ACV  trunks, 
and  equation  [d-19],  times  a geometric  factor  G, can  be  used  tenta- 
tively for  CAB  type  ACV  seal  or  ski  drag.  This  geometry  factor  G is 
less  than  1,  and  depends  on  the  length-to-beam  ratio  of  a CAB  ACV. 

If  we  put  G in  equation  LD-19],  it  takes  a value  of  1 for  the  full- 
peripheral  ACV.  The  factor  G is  tentatively  given  as  1/(1  + i/b). 


EXTERNAL  AERODYNAMIC  DRAG 

The  external  aerodynamic  drag  is  given  by  the  product  of  the  drag 

coefficient  C^  , the  reference  base  cushion  area  S,  and  the  air 
e 

dynamic  pressure  q . The  external  aerodynamic  drag/weight  ratio  is 

2L 

then, 
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RAM  OR  MOMENTUM  DRAG 

The  ram  or  momentum  drag  arises  from  bringing  a constant  mass  flow 
rate  of  air  m from  a velocity  V relative  to  the  ACV  to  a zero  velocity 

relative  to  the  ACV.  Thus, 

®r  ‘ * V ' “ Vj(°c  S«  ) Vo  tD'2' 

where  p V.  is  the  mass  flow  rate  per  unit  area  near  the  ACV  base 

(where  the  flow  is  fully  contracted),  VQ  is  the  vehicle  ve’eoi ty,  and 

A s is  the  fully  contracted  flow  area.  S is  the  daylight  clearance 
-°c  g s 

area,  and  & , the  discharge  coefficient. 

If  the  vehicle  were  fully  cushion  supported,  the  weight  W 
would  equal  wS,  and  the  ram  drag  weight  ratio  would  be 
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If  external  aerodynamic  lift  CL  qg  S provides  partial  ACV  support, 
the  cushion  pressure  p is  less  than  the  specific  weight  w.  Then, 
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DRAG  SUMMARY 

The  total  drag/lift  ratio  D /W  is  given  as  an  appropriate 
sumnation  of  the  previously  derived  component  lift/drag  ratios. 
These  component  ratios  are  repeated  here  for  convenience. 
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1 for  full-peripheral  ACV 
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for  CAB  ACV 
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Foreword 

This  report  is  based  on  a lecture  series  presented  by  the 
authors  at  the  von  Karman  Institute  for  Fluid  Dynamics,  Rhode- 
Saint-Genese,  Belgium,  in  May  1965;  and  at  the  University  of 
Maryland,  College  Park,  Maryland,  in  July  1965.  The  lectures  were 
prepared  under  the  joint  auspices  of  the  David  Taylor  Model  Basin 
and  the  Naval  Air  Development  Center.  They  were  presented  in 
Belgium  under  the  joint  sponsorship  of  the  von  Karman  Institute 
and  the  Advisory  Group  for  Aerospace  Research  and  Development 
(AGARD) ; and  in  Maryland  under  the  sponsorship  of  the  Assistant 
Secretary  of  the  Navy  for  Research  and  Development. 

The  revised  lectures  will  be  presented  as  follows: 

A.  Introductory  Survey 

B.  Air  Cushion  Mechanics 

C.  Internal  Aerodynamics 

D.  Drag 

¥r. — TT  l ii|  C|U  fTTTf  it  inn  for  Ti’rtr'unlf  tTU  ('CAn'l 

F.  Cushion  Contributions  to  Stability 

G.  Seakeeping 

H.  Performance  Summary 


NOTATION 

cushion  length,  ft 
cushion  beam,  ft 

cushion  area,  fta,  measured,  in  plan  view,  to  outer 
edge  of  nozzle  exit 

cushion  perimeter,  ft 

cushion  pressure,  pafg 

total  pressure  of  air  supply  to  cushion  at  nozzle  exit 
or  cushion  entry,  psfg 

flight  velocity  in  earth  axes,  ft/sec 
flight  velocity  in  earth  axes,  knots 
air  density,  slugs/ft3 

daylight  gap  area  (cushion  perimeter  times  mean  daylight 
gap,  h),  ft3 
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nozzle  thickness  parameter  ^ (1  + sin  9)^ 

total  air  volume  flow  rate  through  cushion  system,  ft3 /sec 


cushion  discharge  coefficient 


V p Lp 
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SUMMARY 


A review  is  undertaken  of  some  of  the  more  important  elementary 
relationships  which  determine  the  influence*of  the  air  cushion  on  the 
pitch,  heave,  and  roll  stability  of  ground  effect  machines.  The  analyses 
presented  yield  only  rough  approximations  of  the  quantities  needed  for 
design.  Further,  the  analyses  are  restricted  to  cases  of  vehicles  which 
are  not  in  contact  with  the  surfaces  over  which  they  Crave.)*  whereas 
the  influences  of  surface  contact  are  of  primary  importance  in  most 
practical  problems. 

INTRODUCTION 

The  problem  of  predicting  the  contributions  of  the  cushion  to 
the  stability  characteristics  has  never  been  formulated  in  a fully 
satisfactory  manner.  It  is  regretted  that  it  is  not  practical  to  fully 
rectify  that  unfortunate  situation  at  present. 

The  difficulty  is  twofold.  First,  there  are  some  aspects  of  the 
problem  which  are  not  sufficiently  well  understood  to  support  a fully 
satisfactory  formulation.  In  the  second  place,  even  for  those  aspects 
which  are  reasonably  well  understood,  it  is  found  that  the  different 
elements  of  the  cushion  system  interact  with  one  another  in  a complex 
manner  which  makes  it  very  difficult,  if  not  impossible,  to  arrive  at 
simple  and  reliable  design  formulas.  (There  is  a third  difficulty 
concerning  nonlinearity,  but  it  would  be  idle  to  dwell  on  this  when  we 
cannot  even  deal  properly  with  the  linearized  approximations.) 

The  most  important  contributions  of  the  cushion  to  stability 
involve  the  following  motions: 

Heave  - vertical  translational  motion  (y) . 

Pitch  - rotation  about  the  lateral  axis  (or). 

Roll  - rotation  about  the  longitudinal  axis  (0). 

In  each  of  these  motions,  the  cushion  may  contribute  "stiffness"  (a 
force  or  moment  opposing  the  displacement  from  equilibrium  in  proportion 
to  this  displacement)  and  "damping"  (a  force  or  moment  opposing  the 
motion  in  proportion  to  the  velocity  of  the  motion). 


We  will  ignore  effects  of  deformation  of  the  water  surface  under 
loads  imposed  by  the  cushion.  These  effects  are  important  at  low 
forward  speeds. 

HEAVE  MOTION  WITHOUT  WATER  CONTACT 
The  easiest  of  these  motions  to  understand,  by  far,  is  heave. 

Moreover,  some  of  the  principles  which  govern  heave  motion  are  also  of 
primary  importance  in  the  pitch  and  roll  motions.  We  will  therefore 
consider  the  heave  motion  first,  arid  somewhat  more  carefully  than  the 
other  motions.  Nevertheless,  it  will  be  found  that,  even  in  the  case 
of  heave,  we  will  have  to  content  ourselves  with  a rather  superficial 
treatment  of  the  subject  in  order  to  avoid  an  unduly  laborious  development. 

HEAVE  STIFFNESS 

The  heave  stiffness  coefficient  is  defined  as: 

kh  ' - » <*■> 


This  definition  is  not  mathematically  meaningful,  however,  until  we 
have  specified  which  independent  variables  the  base  pressure  is  con- 
sidered to  depend  upon.  Reflecting  upon  the  relationships  derived  in 
Sections  B and  C,  it  will  be  seen  that  our  problem  would  be  simplest  if 
we  could  write: 

Ap  = Ap(h,  h,  Q) 
or 

Ap  = Ap (h,  h,  pt) 


It  is  easy  to  see,  however,  that  in  heave  motions  of  practical  interest, 

neither  pf  nor  Q are  independent  of  h*  For  example,  assuming  compressor 

rotational  speed  Q constant,  examination  of  Figure  C-4  will  show  that, 

as  h is  decreased  (and  both  f)  S and  \ are  therefore  decreased),  the 

c § 

total  pressure  pt  will  ordinarily  increase,  whereas  Q will  decrease. 
However,  it  may  not  always  be  safe  even  to  assume  Q constant,  because 
the  compressor  torque  will  ordinarily  oscillate  during  heave  oscillation. 
For  a thorough  analysis  one  must  write: 


Ap  = 


Ap(h,  fi,  throttle  setting) 


and  then  derive  a system  of  coupled  differential  equations  describing 
the  behavior  of  engine,  compressor,  and  cushion,  including  an  accounting 
for  the  angular  momentum  of  engine  and  compressor,  and  the  character- 
istics of  the  engine  governor  system,  when  appropriate. 

This  is  not  terribly  difficult  to  do  for  any  one  well-defined 
vehicle,  but  it  is  rather  laborious,  and  the  result  is  so  complex  and 
so  dependent  on  the  peculiarities  of  the  particular  system  considered 
that  it  is  difficult  to  draw  any  simple  conclusions  which  are  of  general 
interest. 

In  the  present  study,  we  will  restrict  ourselves  to  a more  super- 
ficial treatment  of  the  problem  as  follows:  It  was  noted  above  that,  in 

a real  heave  oscillation  (with  no  compressor  stalling),  p^  increases  with 
decreasing  h while  Q decreases  with  decreasing  h;  whereas,  if  h is  held 
artificially  fixed,  Q will  obviously  increase  with  increasing  pfc.  We 
may  surmise  from  this  that  the  correct  value  of  the  stiffness  coefficient 
will  lie  somewhere  between  the  value  calculated  by  holding  Q fixed  and 
the  value  calculated  by  holding  pt  fixed.  We  can  therefore  draw  at  least 
some  qualitative  conclusions  by  comparing  results  of  these  two  calcula- 
tions . 

1.  Q Held  Fixed  - In  this  case  it  is  convenient  to  use  the  following 
relationship  from  Section  B: 


P 2 P ^ he/  c 

This  relationship  is  easy  to  work  with  in  the  two  limiting  cases  x « 1 
and  x » 1.  In  the  first  case, 


and 


4 2 T (1  + sl"  9) 


d(Ap) 


ah 


J Q fixed 


• *<<i 


In  the  second  case, 


and 


2 = Constant 

c 


= - 2 , x » 1 

oh  h 
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These  various  relationships  are  plotted  in  Figure  F-l,  in  the 
form  of  the  nondimens ional  heave-stiffness  parameter  "Sp^S  versus 
nozzle  thickness  parameter  x.  The  qualitative  b 'havior  of  the  "Q-fixed" 
curve  is  indicated  by  a dashed  line.  (This  curve  could  be  determined 
accurately  from  the  relationships  given  in  Section  B,  but  at  the  expense 
of  some  labor;  and,  since  we  are  only  entitled  to  draw  qualitative  con- 
clusions anyway,  it  hardly  seems  worth  while.) 

It  is  seen  in  Figure  F-l  that  the  heave  stiffness  parameter  calcu- 
lated with  Q fixed  covers  a range  from  1.0  to  2.0  as  x increases  from 
zero;  whereas,  with  pt  fixed,  it  covers  a range  from  1.0  to  zero  as  x 
increases  from  zero.  The  most  natural  assumption  to  make  is  that,  with 
throttle  setting  fixed,  the  heave  stiffness  parameter  will  be  relatively 
independent  of  x and  will  have  a value  of  roughly 
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HEAVE  DAMPING 

The  heave  damping  coefficient  is  defined  as 

ch  * • b 4 - s |s  <*■> 

The  same  difficulties  as  to  mathematical  meaning  of  this  difinltion 
(concerning  the  independence  of  variables),  which  were  discussed  above 
for  the  case  of  heave  stiffness,  obviously  apply  here.  As  before,  we 


will  content  ourselves  with  obtaining  rough  estimates  of  the  damping 
I 


(1)  with  Q held  fixed,  and  (2)  with  p£  held  fixed;  and  will  rely  on 

the  fact  that  the  "real"  case  will  lie  somewhere  between  these  ex- 
tremes . 

The  damping  force  arises  primarily  from  the  changing  volume  of 
the  cushion.  The  rate  of  change  of  cushion  volume  is  approximately 

Q'  = -ft  S 

It  is  necessary  to  consider  separately  the  case  ft  < 0 (vehicle  moving 
downward),  when  the  flow  Q'  "squeezes  out"  under  the  normal  efflux  Q; 
and  the  case  h > 0,  when  a portion  Q'  of  the  normal  efflux  Q is  di- 
verted back  into  the  cushion. 


1.  Q Held  Fixed  - We  will  consider  the  four  cases 
and  x » 1»  ft  0 0. 

(a)  x « 1,  h < Q - As  the  flow  Q7  approaches 

r 2 ' 

sure,  it  approaches  the  velocity  / — Ap  and  displaces 


, x « 1,  h > 0 

ambient  pres- 
a thickness 


-5 


The  cushion  pressure  will  be  roughly  the  same  as  if  h were  zero  and 
the  daylight  clearance  were  h - h ; that  is, 


h / h f \ 

Ap  + 6(Ap)  = Ap  ^ J~h'  = & \l  + ~ J 


4 Ap  1 + 


- Ap  hC 
P 


= Ap  1 - 


Q fixed, 
x « 1 » 

h < o 


(b)  x « 1,  ti  > 0 - As  discussed  in  Section  R,  in  the  thin- 
jet  case,  the  momentum  balance  for  the  peripheral  jet  can  be  expressed 


where 


j1  sin  8 + j. 


2 Pt 


is  the  momentum  per  unit  periphery  of  the  jet  leaving  the  nozzle, 
and  j 2 (taken  equal  to  j ^ , in  Section  B)  is  the  net  momentum  per 

unit  periphery  of  the  jet  moving  outward  after  ground  contact.  In 


the  present  case, 


h 4 


Q ~ V "a  Pt  t C 


rrra 


Ap  h = 2 pfc  t 8 In  0+1 


2 h S 


p pt  11  C 


dl ft 


/*£  JT 


a Cap),  ^ 

dh 


/T”-^  fi 


Q fixed, 
, < x « 0, 
h > 0 


(c),(d)  x » 0,  H ^ 0 - These  two  cases  can  be  treated  as 
one.  When  the  jet  Is  very  thick,  ^ i Constant,  Ap  = Pt  » and  the 

base  pressure  can  be  determined  from: 


Q + Q ' = JZ  Ap  hC 


= Q - h S 


. I f Q - h S 
^ = 2 P JT  hr 


2 p > 0©c  hC 


♦ 4, 


Q fixed, 
x » 0, 
h £ 0 
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2.  pt  Held  Fixed 

(a)  x « 1,  h > 0 - In  Section  C,  v/e  saw  that 


Ap  = Pt(2x) 


“Pt  tC 

P 1 


x « 1 


It  therefore  makes  no  difference,  when  x « 1,  whether  one  holds  p^_ 

fixed  or  holds  Q fixed;  the  relationship  between  the  cushion  pressure 
and  h is  the  same.  Therefore,  the  relationships  derived  above  apply 
equally  here;  that  is 


= 

dh 


p fixed, 
x « 1 , 


= 


p fixed , 
x « 1 , 

h > o 


(b)  x « 1,  h ^ 0 - In  this  case  we  have  seen  that  Ap  = p^  • 


Therefore 


'p  fixed, 

m . 0 . J x » 0, 

h $ 0 


An  intuitive  interpretation  of  the  foregoing  limiting-case  rela- 
tionships is  given  in  Figure  F-2.  Note  that,  in  the  case  li  < 0,  the 
situation  is  quite  similar  to  that  found  for  heave  stiffness;  that  is, 
the  range  of  possible  values  of  the  heave  damping  parameter  is  more  or 
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less  symmetrically  disposed  around  the  value  1.0,  and  we  can  be  reason- 
ably confident  that  we  would  ‘not  go  too  far  astray  by  choosing  this 
constant  value  as  an  estimate  for  practical  purposes.  However,  the 
situation  is  considerably  more  complicated  in  tne  case  h > 0.  In  the 
range  x < 1,  the  value  of  the  damping  parameter  is  strongly  dependent 
on  x,  and  can  assume  large  values  when  x is  very  small.  Thus  we  are 
faced,  not  only  with  a situation  in  which  the  damping  parameter  can  vary 
over  a wide  range,  but  also  with  a situation  in  which  (since,  in  general, 
the  damping  coefficient  will  vary  during  a cycle  of  oscillation  as 
passes  through  zero  and  changes  sign)  we  are  not  entitled  to  assume  simple 
harmonic  motion,  even  when  the  amplitude  of  the  motion  is  small. 

Nevertheless,  we  are  faced  with  the  choices  of  (a)  proceeding  on 
the  assumption  that  the  damping  parameter  is  constant,  (b)  undertaking 
a very  complex  investigation  which  takes  into  account  the  variation  in 
damping  parameter,  or  (c)  abandoning  the  investigation  of  heave  motions. 

We  will,  of  course,  follow  course  (a) , assuming,  in  fact,  that  the 
damping  parameter  is  equal  to  1.0;  that  is, 


Note  that,  despite  the  considerable  range  of  uncertainty  involved  in 
this  assumption,  it  is  least  likely  to  give  the  correct  order  of  magnitude, 
under  conditions  of  practical  interest. 

WAVE-EXCITED  HEAVE  MOTION 

Let  us  assume  that  the  vehicle  is  executing  a steady-state  heave 
oscillation,  excited  by  passage  over  a regular  sinusoidal  surface  wave 
system  of  small  amplitude. 

1.  Wave  Length  » Cushion  Length  (\  » l)  - In  this  case  we  can 
measure  the  mean  daylight  clearance  h from  the  bottom  of  the  craft  di- 
rectly below  the  c.g.  to  the  local  water  surface. 

J 
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(Dimensions  distorted  for  clarity) 


The  linearized  equation  of  motion  is: 

u . r 

— y + C,  h + k.  (h  - n)  = 0 

g y h h 

where 

y is  displacement  of  craft  above  mean  position 
h is  distance  of  craft  above  local  water  surface 

h is  distance  of  mean  position  of  craft  above  mean  level  of  water 
surface . 

If  the  wave  system  passes  under^  the  craft  with  frequency  uu 
radians /second , we  can  (with  proper  choice  of  time  origin)  write 

h = h+y-esinuut 
h = y - eui  cos  uut 

The  equation  of  motion  then  becomes 

W 

— y + C.  y + k.  y = k.  e sin  ait  + C,  euo  cos  ujt 
g 7 h ' h h h 
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Now,  assuming  W = Ap  S,  we  can  write  Equations  F-l  and  F-2 


h h 


W S_ 

' hC 


The  equation  of  motion  can  then  be  written 

y + 2 C u>  y + U)  2 y = u)  2 efsin  out  + 2 C—  cos  tut*) 

o o \ 3 0)  / 


where 


r = - / — — 

° 2 V 2 ip  hC 


This  equation  is  easily  solved,  yielding 

f 1 ♦ («ctf  V . r 


i + Ky  + (2  c ^y 

“V,  / V “J 


sin  u)t 


(2ct) 


2 C — 

* UJ 

o 

TTJJ 

1 2 

“o 


Example  - Suppose  h = 0.5  foot;  Ap  = 0.5  l;  -g  = 2,  S = ^jf2 ; and 
C = 3 £ . (These  assumptions  are  roughly  consistent  with  currently 
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operating  vehicles.)  Then  the  motion  is  approximately  as  follows,  for 


craft  of  various  sizes: 
l (ft) 

c 


(rad  tans  /seconds  ) 


25 

0.16 

8.0 


100 

0.32 

8.0 


400 

0.64 

8.0 


Note  the  distinct  resonance  effect,  which  tends  to  become  less  pro- 
nounced for  larger  craft,  and  note  the  tendency  for  the  heave  motion 
to  become  very  small  at  high  frequencies  of  encounter. 

2.  Efffect  of  Wave  Length  - If  the  wave  length  is  not  very  large 
compared  with  the  cushion  length,  it  is  necessary  to  consider  the 
variations  in  daylight  clearance  along  the  length  of  th.e  craft;  that  is, 
instead  of  measuring  h at  the  center  of  the  craft,  we  must  calculate  the 
average  value  of  the  local  gap  h along  the  entire  periphery. 


For  a rectangular  cushion  this  Is  fairly  easy. 


For  a particular  point  along  tne  periphery  we  can  write: 
h^(x,  t)  = h + y - e sin  ^ajt  + 2 n * J 


and  thus  calculate  the  average  gap: 


h(t)  = 


i 1 

b+f.  l 


r2 


h^(x, t)dx 


1 

+ 2 


b 

b+2 


(2  ,t:)  + 2 ‘Ol 


1 

2 


or 


h = h + y 


1/ b 

1+2/b 


sin  (n  j l/\)  1_ 


1 /\ 


cos  (ttA/X)"^  e sin  uut 


(Note:  We  are  averaging  along  the  periphery,  which  is  appropriate  for 

that  part  of  the  excitation  which  arises  from  heave  stiffness.  A smaller 
excitation  term  arises  from  heave  damping,  and  for  this  term  the  average 
should  be  taken  over  the  whole  cushion  area.  We  will  neglect  this  refine 
ment. ) 
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thi3  degenerates  to  the  same  formulation 


Note  that,  in  the  limit 

I 

we  used  previously;  that  is, 

h = ii  + y - c sin  at  , \ » l 

Thus,  our  results  for  the  case  \»  l can  be  used  directly  for  the 
present  case  by  substituting  s , ~^j  for  <t  in  those  results,  where 

the  heave  attenuation  factor  due  to  finite  wave  length  is: 


l/b  sin  (rr  l!  X)  1 

1 1+4/b  tt  l!  \ 1+^/b 


cos  (tt  l/\) 


In  other  words: 


f ampl  (y.)1 

_ 1 

"ampl  (y )' 

€ 

arbitrary  — 

FU  * \) 

e 

[F-4 


The  heave  attenuation  factor  is  plotted  in  Figure  F-3.  For  a craft 
of  practical  length/beam  ratio,  it  is  found  that  there  is  one  range  of 
wave  lengths  slightly  longer  than  the  cushion  length,  and  a second 
range  of  wave  lengths  about  two-thirds  the  cushion  length,  for  which 
the  (leave  motions  will  be  practically  nil  regardless  of  frequency  of 
encounter  uj  . (Indeed,  there  are  an  infinite  number  of  such  ranges 
of  shorter  and  shorter  wave  lengths.  However,  these  further  ranges 
are  of  limited  interest,  because  the  frequency  of  encounter  will 
generally  be  so  high  as  to  preclude  significant  heave  motions  anyway.) 

Before  leaving  the  subject  of  heave  damping,  it  may  be  well  to  * 
review  what  we  have  done. 

1.  8y  making  broad  (and  weakly  justified)  assumptions  regarding 
the  characteristics  of  the  engine-compressor-ducting  system,  we  were 
able  to  obtain  a simple  expression  for  the  heave  stiffness,  Equation 

[f— i 3 • 

2.  With  still  more  sweeping  assumptions,  we  obtained  a simple 
expression  for  heave  damping,  Equation  [F-2]. 
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3.  These  results  enabled  us  (with  the  further  assumptions!  craft 
moving  over  a small-amplitudal  sinusoidal-wave  sea,  traveling  normal 

to  wave  crests)  to  represent  the  heave  motion  by  an  ordinary  differential 
equation  typical  of  a simple  displacement-excited,  damped  spring-mass 
system;  and  to  solve  this  equation,  first  in  the  case  X » l (Equation 
[F-3]);  and  then  in  the  case  of  arbitrary  X (Equation  [F-4]). 

4.  Certain  conclusions  were  indicated: 

(a)  The  natural  frequency  in  heave  is  approximately 


(This  is  borne  out  by  experience.) 

(b)  The  nond imens ional  heave  damping  coefficient  is 


c = 


1 

2 


S_ 

hC 


(It  is  difficult  to  determine  whether  or  not  this  is  in  quantitative 
agreement  with  experience;  but,  at  least  the  effects  of  the  design 
parameters  S/C,  h,  and  Ap  appear  to  be  qualitatively  borne  out  by  ex- 
perience . ) 

(c)  A resonance  effect  is  predicted,  involving  relatively  high 
heave  response  when  the  encounter  frequency  yu  is  nearly  equal  to  natural 
frequency  ujo  , and  a decay  in  heave  response  toward  zero  as  u!  increases 
to  still  higher  values.  (This  is  borne  out  by  experience.  However,  this 
resonance  effect  is  of  little  practical  significance,  because  motions 
under  conditions  which,  according  to  our  analysis,  would  produce  severe 
resonance  effects  are  usually,  in  practical  cases,  dominated  by  water- 
contact  effects.  Resonance  effects  are  seldom  noticeable  in  currently 
operating  vehicles,  for  this  reason.  There  are  also  some  indications 
from  experience  that  "coupling  effects"  between  pitch  motion  and  heave 
motion  can  be  significant.) 

(d)  An  attenuating  effect  depending  on  planform  shape  and  ratio 
of  cushion  length  to  wave  length  is  predicted  (Figure  F-3),  due  to  the 
fact  that  the  craft  tends  to  respond  to  changes  in  average  daylight 
clearance  h,  rather  than  to  changes  at  any  one  point.  (This  is,  at 
least  qualitatively,  borne  out  by  experience.) 
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REMARKS 


The  fact  that  wc  have  obtained  results  which  are  generally  con- 
sistent with  experience  provides  some  encouragement  that  the  main  ideas 
employed  in  the  analysis  were  more  or  less  valid.  It  is  to  be  hoped 
that  the  analysis  itself  has  been  useful  to  the  extent  of  allowing 
some  of  these  ideas  to  be  formalized,  and  some  of  their  consequences 
to  be  explored.  However,  the  reader  is  warned  not  to  expect  anything 
more  from  this  type  of  analysis.  At  almost  every  step  of  the  analysis, 
we  have  resorted  to  simplifying  assumptions  which  far  overstep  the 
bounds  of  what  is  normally  considered  to  be  allowable  and  justifiable. 

It  would  be  a dangerous  mistake  to  suppose  that  the  results  can  be 
safely  applied  for  design  purposes,  or  as  a basis  for  further,  more 
detailed  analyses.  * 

PITCH  AND  ROLL  MOTIONS  WITHOUT  WATER  CONTACT 

The  following  discussion  will  be  pursued  in  as  brief  and  simple 
a manner  as  possible,  attempting  only  to  convey  some  of  the  main  ideas 
involved  in  pitch  and  roll  motions  in  the  absence  of  water  contact, 
without  attempting  to  develop  these  ideas  to  the  point  of  design  appli- 
cation . 

PITCH  MOTIONS 

The  most  common  form  of  attitude  stabilization  is  "compartmenta- 
tion."  The  underlying  idea  is  this:  Suppose  one  took  four  individual 

craft,  each  of  which  was  stable  in  heave,  and  connected  them  together 
with  a rigid  framework.  The  resulting  combination  would  be  stable  in 
pitch  and  roll,  since  each  individual  unit  would  tend  to  seek  and  hold  a 
certain  daylight  clearance  h.  A similar  result  is  achieved  by  introducing 
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PITCH  STIFFNESS  --  Suppose  the  cushion  is  rectangular,  with 
perfect  barriers  which  divide  the  cushion  into  four  equal  compartments 
and  permit  no  flow  whatever  from  one  compartment  to  another.  Then 
each  compartment  is  effectively  an  individual  cushion  with  the  same 


ratio 


S_ 

hC 


(measured  in  level  attitude)  as  the  total  cushion. 


If  the  vehicle  pitches  bow-up  to  a small  anple  a (radians),  the 
mean  daylight  clearance  of  the  fore  and  aft  compartments,  respectively, 
is 


h 


fore 


h + 


= h + 


1 

2 


1 + \ Jt/b 
1 + lib 


l a 


h + 6h 


Applying  our  previous  estimate  for  heave  stiffness  we  obtain 
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■ 


l > 

r • 


> 


6(Ap)fore 
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aft 


The  pitching  moment  is 
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Assuming  W = ApS,  we  obtain 


M_ 

W£ 


i 1 + i ^/b 

8 “I  + l/b 


l 

h 


or 


Now,  to  allow  for  the  fact  that  the  compartmenting  barriers  are  not 

perfect,  let  us  introduce  an  empirical  "barrier  effectiveness  factor" 

E,  and  write 
b 


M_ 

W2 


I E (ll  2 5.^  i 

8 b \ 1 + jJ/b  ) h 


[F-5] 


(It  should  be  pointed  out  that  the  heavV;  stiffness  estimate  introduced 
in  the  derivation  of  Equation  [F-5 ] was  based  on  a qualitative  con- 
sideration of  the  compressor-ducting  characteristics,  which  does  not 
apply  to  the  pitch  and  roll  cases,  unless  each  cushion  compartment  is 
supplied  from  a common  system.  If  all  compartments  are  fed  from  one 
system,  then  the  appropriate  estimate  would  be  based  on  the  "p  Fixed" 
curve  of  Figure  F-l  rather  than  the  "Throttle  Setting  Fixed"  curve,  and 
the  pitch  stiffness  could  be  much  less  than  that  given  by  Equation 
[F-5].) 
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PITCH  DAMPING  — In  the  same  way  that  the  estimate  of  heave  stiff- 
ness was  extended  to  provide  an  estimate  of  pitch  stiffness,  we  can 
extend  the  estimate  of  heave  damping  to  provide  an  estimate  of  pitch 
damping.  However,  that  exercise  will  be  omitted,  since  the  principles 
involved  should  be  obvious,  by  this  time,  to  the  serious  reader;  and 
the  result  is  of  rather  limited  practical  value  for  the  various  reasons 
discussed  earlier. 

PITCH  NATURAL  FREQUENCY  --  We  will  consider  pitch  motion  only  to  the 
extent  of  estimating  the  natural  frequency.  This  we  can  do  by  con- 
sidering the  equation  of  motion  for  free  oscillation  in  pitch,  neglecting 
damping: 


, ..  m 

I a ■ t 

y da 


a = 0 


From  Equation  [F-5 ] we  have 


i „ + * 


M ± - I E 1—1—  w 1 

da  ~ 8 cb  1+  i/b  "h 


Also,  we  can  write  the  moment  of  inertia  about  the  lateral  axis  in  the 


T W - 2 

I = - r 

y g y 

where  r is  the  "radius  of  gyration"  about  the  lateral  axis.  The 

y 

equation  of  motion  becomes 


s + 1 E a + 1 —(a.)  s a 

“ 8 Eb  v 1 + i/b  ) \r  / h 


and  the  solution  is: 


1 U)  t 
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where 


UJ 

o,a 


fEb  1 + j i/b 
B~  ~1  + i/b 


(Example  - Suppose 


JL 

2 ’ 


i/b  ■=  2, 


r 


y 


i 
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then 


uu  = 0.61  = 0.61  uj 

o, or  V h o 


0-6] 


This  result  is  roughly  consistent  with  experience.  Most  experiments 
have  indicated  the  natural  frequency  in  pitch  to  be  between  half  to 
three-quarters  of  the  natural  frequency  in  heave,  .) 


ROLL  MOTIONS 

ROLL  STIFFNESS  --  The  development  of  an  expression  for  roll  stiff- 
ness follows  exactly  the  above  development  for  pitch  stiffness,  with 
the  roles  of  length  and  beam  interchanged.  The  result  can  be  v'ritten 
down  by  inspection  from  Equation  [F-5],  substituting  rolling  moment  £ 
for  pitching  moment  M;  roll  angle  0 for  pitch  angle  a>  beam  b for 
length  i,  and  vice  versa. 


Z_ 

Wb 


Eb 

8 \ £/b  + 1 


b 
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[F-7] 


ROLL  NATURAL  FREQUENCY  Similarly,  the  result  for  natural  frequency 

in  roll,  a)  , can  be  written  down  by  inspection  from  Equation  [F-6]. 

o,0 


U) 

0,0 


i/b  + ^\ 
i/b  + 1 


b_ 

r 
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[F-8] 


ATTITUDE  STABILITY  WITHOUT  C OMPA RTME NTA T I ON 

The  modern  design  trend  for  air  cushion  vehicles  is  toward  ex- 
tremely small  daylight  clearance,  accepting  some  degree  of  water 
contact  at  practically  all  times.  Under  these  conditions,  it  becomes 
comparatively  unimportant  precisely  what  degrees  of  pitch  and  roll 
stiffness  the  cushion  provides  in  the  absence  of  water  contact.  A 
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serious  question  arises  as  to  whether  the  design  complications  and  other 
penalties  associated  with  compartment ing  jets  or  baffles  are  really 
j ustif led • 

On  the  other  hand,  it  is  certainly  most  undesirable  to  have  the 
craft  become  actively  unstable  in  pitch  and  roll,  in  those  rare  in- 
stances when  the  craft  might  operate  free  of  contact  with  the  surface. 

The  craft  might  become  practically  unmanageable  under  such  circums --antes 
unless  positive,  fast-acting  pitch  and  roll  controls  were  provided  - which 
in  themselves  are  usually  unwanted  design  complications. 

Peripheral  jet  type  craft  are  ordinarily  actively  unstable  unless 
fitted  with  compartmenting  devices.  This  was  especially  true  of  the 
old-style  flat-bottom  designs,  which  were  subject  to  violent  destabi- 
lizing effe.cts  of  the  "cross  flow"  of  air  through  the  cushion  from  the 
low-side  jet  toward  the  high  side.  Under  certain  conditions,  a "venturi 
throat"  type  of  flow  could  establish  itself  along  the  low  side  of  the 
base,  forming  a low-pressure  region  which  sucked  the  low  side  down  to  the 
surface. 


The  modern  peripheral- jet  design  trend  toward  deep  flexible  trunks, 
with  deeply  recessed  base  in  between,  should  presumably  greatly  alleviate 
these  cross- flow  effects,  but  apparently  compartmentation  is  still  re- 
quired to  avoid  some  degree  of  attitude  instability. 

Plenum  type  craft,  however,  usually  exhibit  a reasonably  acceptable 
degree  of  pitch  and  roll  stiffness,  even  if  compartmenting  devices  are 
omitted  from  the  design.  There  are  several  phenomena  which  can  contribute 
to  attitude  stability  of  plenum  type  craft,  but  we  will  consider  only 
the  most  important  one,  which  involves  pressure  forces  acting  near  the 
lower  edge  of  the  plenum  wall. 


***>*■» 


The  efflux  through  the  orifice  of  height  h "necks  down"  to  a vena 
contracta  of  height  h . It  is  easy  to  show  that  there  must  be  a 
reaction  force  per  unit  length  of  the  plenum  of  2 Ap  h on  the 
whole  plenum  system.  A portion  Ap  h of  this  reaction  is  readily 
identified  as  an  unbalanced  pressure  reaction  on  the  far  vertical 
wall  of  the  plenum  opposite  the  orifice.  The  remainder  is  ascribed 
to  a "lip  force"  F^  , identified  with  the  pressure  distribution  near 
the  orifice  lip.  (Note  that  the  lip  force  F is  a hypothetical  force 

AJ 

in  the  sense  that  the  pressures  on  the  lip  actually  act  in  the  opposite 
direction  from  F.  . The  Force  F.  is  the  vector  difference  between  the 

l l 

actual  lip  force  and  the  lip  force  which  would  act  if  the  pressure  were 
uniformly  equal  to  Ap.)  The  magnitude  of  F^  is  determined  from 


-p  cos  9 + Ap  h 


2 Ap  ft  h 
c 


It  , , „ - 1 

/ 15  cos  9 
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where  l , in  the  present  context,  is  the  length  (normal  to  the  plane 
of  the  sketch)  of  the  plenum. 

Using  the  relation  from  Section  B: 


J)  i 1 i + F.?-9  0 I 

c ^ "-+'  \ (1  + sin  0)  - sin  0 cos  0 j 

we  get 


n + 2 
TT  “ 2 


(1  + sin  0)  - sin  0 cos  & 


We  are  now  in  position  to  calculate  the  roll  stiffness  of  a 
plenum  type  vehicle  with  t/b  » 1. 


c.g.  Height  c 


Referring  to  the  sketch,  we  have 

£ = - (f  . , - F „ G cos  0 + F i " F . - ! si 

\ Jl, l i, 2 i \ j?,1  1,2  2 


,1  " *Plh  + *2 


n - 2 


(1  + sin  9)  - sin  0 cos  0 


Fi, 2 4 lp(h  • 


(1  + sin  0)  - sin  0 cos  0 


F - F 
1,2  l, l 


AP  0 b l 

n-  + (1  + sin  9)  - sin  0 cos  0 


or,  if  W = Ap  b jt 


£_ 

Wb 


P 1 

— cos  0 - — sin  0 

— 0 , 1/ b » 1 [F-9] 

(1  + sin  0)  - sin  9 cos  9 


This  relationship  is  illustrated  in  Figure  F-4.  According  to  this 
formula,  the  plenum  type  craft  will  always  have  stable  roll  stiffness 
when  9 is  zero  or  negative,  becoming  less  stable  with  increasing  0 until 
it  becomes  actively  unstable  when  9 is  so  large  that  the  vectors  F^ 
point  above  the  c.g. 

Note  that  roll  stability  (without  water  contact)  improves  with  in- 
creasing c.g.  height  or  decreasing  beam,  a slightly  surprising  result, 
perhaps.  (However,  a high  ratio  of  c.g.  height  to  beam  G/b  is  very 
dangerous  in  practice,  when  water  contact  effects  are  present.) 

When  vehicles  of  practical  length/beam  ratios  are  considered,  it 
is  found  that  Equation  [F-9]  is  approximately  valid  if  9=0.  The  craft 
will  be  less  stable  than  Indicated  by  Equation  [F-9]  if  9 > 0,  more 
s table  if  0 <0. 

When  pitch  is  considered,  a result  for  the  case  b /£  « 1 can  be 

written  down  by  inspection  from  Equation  [F-9]  (replacing  £ with  M and 

» 

b with  l) ; but  this  case  is  of  little  practical  interest.  However,  we 
can  write  down  a useful  formula  for  the  isolated  case  9=0,  in  which 
case 


M_  i tt  - 2 G 

W l ~ ‘ TT  + 2 l 0 


9=0° 


[F-101 


The  above  developments  have  been  based  on  the  assumption  chat  the 
walls  of  the  plenum  are  sufficiently  deep  compared  to  the  gap  h that 
the  pressure  drop  occurring  near  the  lip  is  felt  almost  entirely  on 
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I 

I 


I 


i 


the  walls,  and  not  on  the  base  of  the  vehicle.  If  the  walls  are 
shallow  compared  to  the  gap  h,  so  that  the  pressure  drop  Is  felt 
on  the  base,  then  the  vehicle  can  be  much  more  stable;  but  of  course 
there  Is  a corresponding  loss  of  efficiency  of  the  plenum  In  pro- 
vidirtg  lift. 

We  saw  in  Section  B that  there  is  no  clear  line  of  demarcation 
between  the  aerodynamics  of  peripheral  jet  craft  and  plenum  craft.  The 
latter  can  be  thought  of  as  simply  a form  of  the  former,  but  with  very 
thick  jets. 

By  the  same  token,  it  must  be  realized  that  the  same  consider- 
ations which  led  to  Equations  [F-9]  and  [F-10]  for  plenum  craft,  apply 
also  (in  modified  form)  to  ordinary  peripheral- jet  craft,  becoming 

more  and  more  important  as  the  nozzle  thickness  parameter,  x = --(l+sinQ), 
is  increased. 

REMARKS 

The  remarks  made  at  the  end  of  the  study  of  heave  motions  apply 
equally  here;  that  is,  the  derivations  of  formulas  concerning  pitch  and 
roll  motions  were  intended  primarily  as  a means  of  formalizing  some 
main  ideas,  and  exploring  some  of  their  consequences.  The  formulas  are 
not  likely  to  be  reliable  for  design  application. 

It  is  even  more  important  not  to  lose  sight  of  the  following 
fact:  The  whole  of  this  section  has  concerned  itself  with  motions  in 

the  absence  of  water  contact.  Modern  design  trends  seem  to  point 
toward  practically  all  types  of  GEM's  operating  in  intimate  contact 
with  the  water  (or  whatever  surface  they  are  traveling  over)  at  practi- 
cally all  times.  The  contributions  of  the  cushion  to  stability  ’(except, 
perhaps,  to  heave  stability)  thus  become  much  less  significant  than 
the  forces  and  moments  imposed  by  contact  between  the  water  (or  other 
surface)  and  the  trunks,  sidewalls,  seals,  hull,  or  whatever  portions 
of  the  vehicle  experience  the  contact. 


-25 


I 


These  extremely  important  effects  are  discussed  (not  adequately, 
by  any  means,  but  to  the  best  of  the  authors'  ability)  in  a later 
section. 


Aerodynamics  Laboratory 
David  Taylor  Model  Basin 
Washington,  D.  C. 
February  1967 
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Figure  F-4  - Roll  Stiffness  for  a Plenum  Type  Craft 
Without  Compartmentation 
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Foreword 

This  report  is  based  on  a lecture  scries  presented  by  the 
authors  at  the  von  Karman  Institute  £or  Fluid  Dynamics,  Rhode- 
Saint-Cenese , Belgium,  in  May  1965;  and  at  the  University  of 
Maryland,  College  Park,  Maryland,  in  July  1965.  The  lectures  were 
prepared  under  the  joint  auspices  of  the  David  Taylor  Model  Basin 
and  the  Naval  Air  Development  Center.  They  were  presented  in 
Belgium  under  the  joint  sponsorship  of  the  von  Karman  Institute 
and  the  Advisory  Group  for  Aerospace  Research  and  Development 
(AGARD) ; and  in  Maryland  under  the  sponsorship  of  the  Assistant 
Secretary  of  the  Navy  for  Research  and  Development. 

The  revised  lectures  will  be  presented  as  follows: 

A.  Introductory  Survey 

B.  Air  Cushion  Mechanics 

C.  Internal  Aerodynamics 

D.  Drag 

— Beau  Dpt  imizAtion  f ot^&i  d-awall — GKM_ 

F.  Cushion  Contributions  to  Stability 

G.  Seakeeping 

H.  Performance  Summary 
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NOTATION 


local  wave  elevation,  feet 

amplitude  of  sinusoidal  wave  component,  feet 

stationary  frequency  of  sinusoidal  wave  component, 
radians/sec  (frequency  with  which  wave  propagates 
past  a fixed  point) 

frequency  of  encounter  of  sinusoidal  wave  component, 
radians/sec  (frequency  with  which  wave  is  encountered 
by  a moving  point) 

wave  length,  feet 

distance  along  horizontal  axis  in  direction  opposite 
to  wave  propagation,  feet 

acceleration  of  gravity,  ft/sec2 

velocity  of  wave  propagation,  ft/sec 

mass  density  of  water,  slugs/ft3 


arbitrary  integer 
mean-square  of  h,  ft* 


V wind  velocity,  knots 

w 

[A(y.)]2  spectral  density  of  wave  system,  as  function  of 

stationary  frequency  of  components,  ft3-sec 

[A(uj)]2  spectral  density  of  wave  system  as  function  of  frequency 

of  encounter  of  components,  fts-sec 

H wave  height  of  composite  wave,  feet 

H mean  of  all  values  of  H 

av 

mean  of  the  highest  one-third  values  of  H 
mean  of  the  highest  one-tenth  values  of  H 

V velocity  of  moving  point  (craft),  ft/sec;  taken  positive 

in  direction  opposite  to  wave  propagation 

a pitch  angle,  radians 


heave  displacement,  feet 


lit 


i 


m 

oar 


pitch  magnification  factor  for  linear  response  to  a 
sinusoidal  wave 

heave  magnification  factor  for  linear  response  to  a 
sinusoidal  wave 

2 

mean-square  value  of  or,  radians 

mean-square  value  of  y,  ft2 

cushion  area,  ft2 

cushion  length,  feet 

cushion  beam,  feet 

mean  cushion  length  (S/b),  feet 
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SUMMARY 


The  analytical  and  experimental  tools  available  for  application 
to  GEM  seakeeping  problems  are  reviewed  briefly.  The  one  available 
set  of  seakeeping  data  for  a full-scale  GEM  is  reviewed,  and  is  found 
to  constitute  encouraging  evidence  that  GEM  seakeeping  problems  will 
be  tractable. 

INTRODUCTION 

The  term  "seakeeping"  refers  to  the  riding  qualities  of  a vehicle 
which  travels  on  a "seaway";  that  is,  in  the  ocean-surface  environment. 

A vehicle  with  good  seakeeping  qualities  moves  through  or  over  the 
ocean  waves  with  a minimum  of  the  violent  motions  which  would  cause  dis- 
comfort to  the  passengers  and  interfere  with  the  functions  of  crew  and 
equipment,  a minimum  of  the  structural  loads  which  require  additional 
investment  in  structural  weight  at  the  cost  of  payload,  and  a minimum 
of  the  additional  resistance  which  causes  degradation  of  speed  and  range. 

The  question  of  the  seakeeping  characteristics  of  Ground  Effect 
Machines  is  probably  the  most  crucial  of  the  questions  which  will  affect 
the  future  economic  and/or  military  importance  of  these  vehicles.  It  is 
not  a simple  question,  of  course,  and  one  cannot  expect  a simple  answer. 
Even  for  displacement  ships,  the  progress  toward  scientific  methods  of 
predicting  the  seakeeping  qualities  of  a given  hull  form  has  been  com- 
paratively recent;  and  much  still  remains  to  be  done  in  the  problem  of 
taking  proper  account  of  seakeeping  in  the  initial  design.  The  science  of 
ship  design  is  one  of  the  oldest  sciences,  yet  still  depends,  in  many 
aspects,  on  gradual  evolution  and  the  "designer's  eye." 

The  very  short  history  of  Ground  Effect  Machines  has  permitted 
only  the  beginnings  of  an  evolutionary  process,  and  only  a very  few 
designers  (mostly  in  the  British  Hovercraft  projects)  can  lay  claim  to 
even  the  most  rudimentary  bases  for  development  of  a "designer's  eye." 
Nevertheless,  this  short  history  has  already  produced  encouraging  evidence 
that  the  seakeeping  problems  of  GEM's  will  be  found  tractable,  to  a 
degree  that  few  would  have  dared  to  hope,  before.  This  early  evidence 
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of  tractability  pertains  largely  to  the  remarkable  cushioning  effects 
of  the  air  cushion  itself,  plus  the  remarkable  early  successes  in  the 
development  of  flexible  understructure. 

Unfortunately,  the  important  information  concerning  seakeeping 

/ 

of  GKM's  has  usually  found  its  way  into  the  literature  in  very  frag- 
mentary forms.  Despite  the  predominant  importance  of  this  subject,  the 
following  discussion  will  have  to  be  limited  to  a very  brief  review  of 
the  most  significant  of  the  available  mathematical  methods,  plus  a rather 
superficial  examination  of  the  most  readily  accessible  experimental 
evidence. 

DESCRIPTION  OF  NATURAL  WAVES 

The  following  material  is  intended  as  an  introduction  to  the 
simpler  aspects  of  this  complex  subject,  for  benefit  of  those  readers 
who  are  unfamiliar  with  the  problems  of  mathematical  description  of 
the  sea. 

Unless  otherwise  stated,  this  material  is  based  on  Vossers1  "Fun- 
damentals of  the  Behavior  of  Ships  in  Waves"  (Reference  1). 

SIMPLE  WAVES 

The  simplest  non-trivial  solution  for  motion  of  the  surface  of  a 
body  of  water  is  the  one-dimensional  sinusoidal  wave  propagating  in  the 
(negative,  by  choice)  direction  of  the  x-axis.  The  wave  elevation 
measured  from  the  mean  is 


h = h s in  i ^ x + p,t^ 

where  R is  the  amplitude  of  the  wave  (half  the  wave  height  H) , X is  the 

t 

wavelength,  and  is  the  frequency  (in  radians/sec)  with  which  the 
wfive  propagates  past  a fixed  point. 


i 


If  the  wave  "steepness"  Is  sufficiently  snail  ^ S "fo)  * an<*  t*'e 
water  sufficiently  deep  (depth  2:  , the  frequency  is  given  by 


^ = V2 tt  g/X  , radians/sec 
The  velocity  of  propagation  is  obtained  from 


£ 

X 


or 


c = 


1L 

2tt 


'££  = £ 
2n  y, 


and  the  energy  stored  in  the  wave  motion,  per  unit  area  of  the  sea 
surface,  is 


E = -r  p h 
2 Kw 


(This  is,  in  other  words,  the  energy  which  would  be  required  to  estab- 
lish the  wave  motion,  by  a conservative  process.) 

Natural  waves  are  not  precisely  sinusoidal.  The  trochoidal  wave 
form  is  used  for  some  purposes,  as  a somewhat  better  approximation  of 
natural  wave  form.  However,  the  deviations  from  sinusoidal  are  very 
slight,  except  in  cases  of  extremely  steep  waves,  and  the  sinusoidal 
representation  is  far  more  convenient  for  purposes  of  mathematical 
analysis. 

IRREGULAR  WAVES 

The  sea  displays  many  characters.  Some  of  the  common  terminologies 

are: 

Wind  Sea  — A sea  which,  formerly  calm,  is  building  up  under  the 
action  of  winds;  characterized  by  a high  degree  of  irregularity. 

Swell  — A sea  in  process  of  subsiding,  in  relatively  calm  wind; 
comparatively  regular,  often  comparatively  long-crested. 

Fully  Arisen  Sea  — A sea  which  has  been  exposed  to  a more  or  less 
constant  wind  condition  for  some  time,  and  has  reached  a condition  of 
equilibrium  between  addition  of  energy  from  the  wind  and  dissipation  of 
energy  from  the  viscous  action  of  the  water. 
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Long-Created  Sea  — Waves  all  propagating  In  generally  the  same 
direction,  so  that  relatively  long,  parallel,  unbroken  crests  are 
visible.  (Sometimes  used  to  denote  mathematically  one-dlmenslonal 
waves,  which,  of  course,  do  not  occur  from  natural  causes.) 

Short-Crested  Sea  (Confused  Sea)  — Waves  propagating  In  a variety 
of  directions,  so  that  crests  are  short,  broken,  and  distinctly  three- 
dimensional  In  nature. 

The  most  fundamental  property  of  the  sea's  surface  is  randomness. 

r . 

In  principle,  It  would  be  possible  to  mathematically  reproduce  the 
polnt-by-point  motion  of  any  finite  area  of  ocean  surface  over  any 
finite  span  of  time.  The  futility  of  doing  so  is  apparent  upon  re- 
flecting that  this  particular  motion  will  not  occur  again,  anywhere, 
in  a million  years.  However,  certain  statistical  properties  of  the  sea 
(such  as  the  mean  wave  height,  the  relationship  between  the  mean  height 
and  mean  frequency,  etc.)  are  repeated;  and  it  is  the  statistical  prop- 
erties which  we  must  seek  to  represent  and  understand. 

The  statistical  properties  of  virtually  any  sea  condition  can  be 
satisfactorily  represented  by  superposition  of  an  infinite  number  of 
simple  sinusoidal  waves,  each  of  infinitesimal  amplitude,  propagating 
in  a certain  random  distribution  of  directions  with  random  phase. 

For  one -dimensional  (uni-directional)  waves,  this  representation 
can  be  expressed 


h 


N 

Jj  Ri sln  [17 x f “1 ' + ei] 


where  N is  an  arbitrarily  large  number,  and  where  the  phase  angle  is 
random,  with  uniform  distribution  over  the  interval  0 < e(u)  < 2n  . 

(It  might  seem  natural  at  first  glance  to  let  N -*  00  and  replace  the 
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summation  by  an  Integral;  bu-  if  this  Is  attempted  by  elementary 
methods,  the  resulting  integral  is  found  to  be  indeterminate.)  Since 
each  component,  owing  to  the  random  phase,  is  a random  quantity, 
the  sum  h,  for  any  given  value  of  x and  t,  is  also  a random  variable, 
of  .mean  value  zero;  and,  according  to  the  central  limit  theorem  of 
probability  theory,  has  a normal  (Gaussian)  distribution.  That  is, 
for  arbitrary  choice  of  x and  t,  the  probability  that  h will  lie  L^cween 
values  a and  b is 


P (a  < ii  < b) 


1 

,/  2n  m 
V o 


f'2"o 


■i? 


i 


j 


| 


E 


a b 


where  m is  the  variance  of  the  random  variable  h;  or,  since  the  wave 
o 

elevation  h has  a moan  value  of  zero,  m^  is  the  mean-square  value  of 
h.  Since  the  mean-square  amplitude  of  the  sum  of  a series  of  sine 
waves  of  different  frequencies  is  half  the  sum  of  squares  of  the  com- 
ponent amplitudes,  we  have 


N 

1 V 
21^ 


h 


2 

i 


-s- 


m ■ 
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Now,  there  la  no  difficulty  In  converting  the  summation  of  aquared 
amplitudes  to  Integral  form.  Lettering  N -*  • , we  define 


a v - 2 

[A(p,)]  dp,  ■ l_ 

dp 

where  the  summation  is  carried  out  for^Ll  compontnts  h^  In  the  fre- 
quency interval  from  p to  p + dp  . Then 


®0  “ \ | [A(p)]  dp 

o 

Evidently,  a great  deal  of  Information  is  available,  concerning 

2 

the  character  of  the  sea,  if  the  function  [A(p)]  is  known.  This 
function  is  called  the  "energy  spectrum"  (owing  to  the  fact  that  the 
energy  associated  with  a simple  waye  is  proportional  to  the  square  of 
its  amplitude),  and  the  value  of  the  function  for  a given  frequency 
is  called  the  "spectral  density."  The  important  science  of  forecasting 
ocean  waves  is  based  on  empirical  knowledge  of  the  dependence  of  the 
energy  spectrum  on  the  intensity  and  duration  of  winds  and  on  the  "fetch" 
(length  of  ocean  surface,  in  the  upwind  direction,  over  which  the  wind 
prevails).  ' 

The  best-known  of  the  mathematical  representations  of  ocean 
energy  spectra  are  those  of  the  "Neumann  sea,"  representing  fully  arisen 
seas  corresponding  to  arbitrary  wind  velocities  , in  knots 
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These  Neumann  spectra  are  semi-empirical  representations  of  time- 
histories  of  the  wave  action  at  a fixed  point  in  a real  fully  arisen 
sea,  which  is,  of  course,  not  one-dimensional.  They  are  frequently 
used  to  represent  one-dimensional  seas,  however,  which  in  turn  are 
assumed  to  be  representative  (after  proper  adjustments  for  the  speed) 
of  the  sea  "sensed"  by  a craft  moving  with  or  against  the  dominant 
direction  of  wave  propagation. 

With  a knowledge  of  the  energy  spectrum,  and  knowledge  that  the 
distribution  of  wave  elevations  is  Gaussian,  numerous  statistical 
properties  can  be  calculated.  For  the  Neumann  spectra,  for  example, 
it  is  found  that: 

a.  Average  wave  height,  measured  between  successive  crest  and 
trough. 


(Since  J m la  the  root-mean-square  value  of  the  wave  elevation  h, 
o 

one  would  correctly  guess  that  this  should  be  about  half  the  average 
crest-to-trough  dimension  H .) 

HV 
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b.  "Significant  wave  height,"  defined  as  the  average  of  the 
highest  one-third  of  the  wave  heights, 

hi/3  ■ 4y^o 

c.  "Maximum"  wave  height,  defined  as  the  average  of  the  highest 
one-tenth  of  the  wave  heights, 

Hl/10  " 5ml^  mo 

These  and  other  statistical  properties  of  fully  arisen  seas  are 
given  In  a very  useful  tabulation  compiled  by  Wilbur  Marks,  and  repro- 
duced here  as  Table  G-l. 

ENERGY  SPECTRA  AS  FUNCTIONS  OF  FREQUENCY  OF  ENCOUNTER 

The  energy  spectra  were  considered  above  as  functions  of  the  circular 
frequencies  p,  with  which  the  simple  wave  components  propagate  past  a 
fixed  point. 

A simple  wave  of  frequency  p,  propagates  with  velocity  c = & . 
Therefore,  a point  moving  normal  to  the  crests  at  velocity 
encounters  this  wave  at  frequency 
« 

c + V , V v 

s— 2 E + V‘Y'1 

Where  V la  considered  positive  for  a point  moving  in  the  opposite 
o 

direction  to  wave  propagation  (head  sea)  and  negative  for  a point  moving 
In  the  same  direction  as  wave  propagation  (following  sea). 

Now  wa  can  express  the  energy  spectrum  [A(p)]8  in  terms  of  frequency 
of  encounter  uj,  taking  care  that  the  "energy"  contained  in  the  frequency 
interval  u>  to  uj  + duu  Is  the  same  as  that  In  the  corresponding  interval 
p to  p + dp  . 

[A(p)]8  ■ LA  (a,)]8  du> 

■ [A(uu)]8  (l  + 2p  ~)  dp 
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[A  ( cu)  ] 


U(p)f 


1 + 2p.  Vo/>g 


H (1  + M.  VQ/g) 


In  the  case  of  a head  sea  V > 0 , the  energy  is  spread  out  over  a 

o 

wider  band  of  higher  frequencies  than  for  the  s tationary-point  spectrum; 
whereas  for  following  sea,  at  moderate  negative  velocities,  energy  tends 
to  be  concentrated  in  narrower  bands  at  low  frequency,  the  spectral 
density  becoming  infinite  at  the  frequency  of  encounter 

» - -15- 

max  4 V 

o 

(Note:  V0  < 0,  so  u>  > 0).  This  singularity  in  the  spectral  density 

is  integrable,  of  course,  so  that  the  energy  in  any  given  band  of  fre- 
quencies is  finite.  The  following-sea  spectra  are  double-valued  in  the 
range  of  frequencies  for  which  it  is  possible  to  encounter  waves  of  two 
wavelengths  at  the  same  frequency  of  encounter,  and  extend  into  the 
"negative  frequency"  range.  (The  physical  interpretation  of  "negative 
frequency"  for  the  following  sea  is  merely  that  the  waves  are  propa- 
gating with  velocity  c < |Vq|  , and  are  being  overtaken  by  the  moving 
point.  The  negative  spectral  density  is  seen  to  pose  no  problem  when 
it  is  considered  that,  in  integrating  the  energy  content  in  any  frequency 
band,  a path  of  increasing  values  of  p is  followed,  so  that  the  product 
[A(uj)]S  du;  is  always  positive.  For  example,  one  can  write 


1 * 

2 J [A(n)  ^dp 


max  -oo 

“ 2 J tA(u))]2  du)  + J [A(uj)]2  dcu 

O 0) 


In  the  final  form  of  the  equation,  the  first  integration  is  carried 

out  from  left  to  right  In  the  upper  half  of  the  plane  ^ [A(uu)  ]2^ versus  uj  ; 
and  the  second  is  carried  out  from  right  to  left  in  the  lower  half 
plane.) 
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Spectral  Density  Versus  Frequency  of  Encounter 
(State  3 Sea.  One-Dlmenslonal  Idealization) 


Some  examples  are  given  in  the  above  sketch  for  a state  3 sea  (see 
Table  G-l)  and  for  speeds  appropriate  to  slow  surface  vessels.  It  should 
be  pointed  out  that  when  one  goes  to  much  higher  speeds,  there  is  little 
practical  difference  between  head  and  following  sea.  The  respective 
frequency  spectra  resemble  reflections,  about  the  origin,  of  one  another. 
The  following-sea  spectrum  still  has  the  singular  "spikes”  near  the  origin, 
of  course;  but  they  contain  practically  no  energy. 
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PREDICTION  OF  MOTIONS  BY  LINEAR  SUPERPOSITION 
Consider  the  motions  of  a surface  craft  traveling  normal  to  the  wave 
crests  of  a simple  sinusoidal  wave  system.  Ordinarily,  we  expect  oscil- 
lations in  pitch  (a)  and  heave  (y)  to  be  induced  by  passage  over  the  wave. 

If,  for  arbitrary  wave  length  and  frequency  of  encounter,  (1)  the 
oscillations  In  pitch  and  heave  are  sinusoidal,  and  (2)  the  amplitudes  of 
oscillation  are  directly  proportional  to  the  amplitude  of  the  wave,  then 
it  is  possible  to  represent  the  response  to  an  irregular  wave  as  the  sum 
of  the  responses  to  an  arbitrarily  large  number  of  simple  waves  of  which 
the  irregular  wave  is  assumed  to  be  composed. 

In  other  words,  if  the  elevation  of  the  simple  wave  is  expressed: 


h ” h sin  u/t 

then,  in  order  for  linear  superposition  to  be  valid,  the  motions  must  be 
expressible  as  follows: 


o <1 

y 


M sin  /"out  + 0 ^ 

My  K 8ln  (*t  + 0y) 


where  the  "magnification  factors"  and  My  and  relative  phase  angles  0_^ 
and  0y  may  depend  on  any  or  all  of  the  variables 

V craft  velocity 

o 

\ wavelength 

u>  frequency  of  encounter 

but  must  not  depend  on  the  wave  amplitude  h.  (The  pitch  and  heave  motions 
need  not  be  independent;  that  is,  heave  displacement  may  induce  pitching 
moment,  etc.,  so  long  as  the  above-mentioned  conditions  are  met.) 

In  this  case,  if  we  denote  the  mean-square  amplitudes  of  wave  eleva- 
tion, pitch  angle,  and  heave  displacement  as:  m^  (feet2),  m (radian2), 


and  m (feet3),  then  for  a given  craft  moving  at  given  velocity  VQ  over 
an  Irregular  wave  composed  of  N simple-wave  components,  we  have 


if  t- 
1 


- 7? 


■ 11(4 


Note  that,  for  given  velocity  Vq,  the  magnification  factors  can  be  con- 
sidered to  be  functions  of  ou  only  (even  though  they  may  depend  physically 
upon  X),  since  (given  Vo  and  id)  , X Is  determined.  It  Is  determined  uniquely, 
In  fact, • provided  one  properly  Identifies  which  of  the  two  possible  values 
of  X applies  in  cases  of  following  sea*  positive  id.  (Refer  to  the  pre- 
ceding discussion  and  sketch,  concerning  expression  of  the  energy  spec- 
trum of  the  sea  as  a function  of  frequency  of  encounter  (page  10). 

The  sums  can  be  converted  to  integrals,  letting  N — For  head  sea: 


J [A (id)  ]"  duu 


7T I [M*H  [A(u,)]  duj 


co 

r 

moy  ” 2 (*•*)]  [a(U))J  d< 
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(We  have  already  considered  the  expression  for  ra  , in  the  case  of  a 

o 

following  sea,  as  the  sura  of  two  Integrals.  The  expressions  for  m 


OQf 


and  rn^  are  directly  analogous.  However,  we  will  confine  the  following 


discussion  to  the  case  of  a head  sea.) 
The  functions 


{[■„(  u>)j  [a(w)]  | and  |[My(w)]  [a  ( oj)  J 


V are  called 


the  energy  spectra  of  the  pitch  and  heave  motions.  They  are  obtained  by 
the  simple  expedient  of  multiplying  the  energy  spectrum  of  the  Sc«  by 
the  squares  of  the  respective  magnification  factors.  From  them,  it  is 
possible  to  determine  innumerable  useful  statistical  properties  of  the 
motions  (average  amplitude  of  oscillation,  average  of  the  one-third 
highest  oscillations,  probability  of  exceeding  a certain  displacement  in 
a certain  time  period,  etc.).  (Note  that,  regardless  of  what  phase 
relationships  may  exist  between  the  elementary  sinusoidal  components  of 
the  wave  and  the  corresponding  elementary  sinusoidal  components  of  the 
motions,  if  the  wave  components  have  random  phase,  then  the  motion  com- 
ponents will  have  random  phase.  Therefore,  the  composite  motion  displace- 
ments will  have  Gaussian  distribution,  with  variance  m and  m , re- 
oar oy 

spectlvely . ) 

Two  problems  remain: 


1.  Justification  of  the  assumption  that  the  responses  to  simple  waves 
are  linear  (i.e.,  sinusoidal  with  amplitudes  proportional  to  wave  ampli- 
tude). 

2.  Determination  of  the  magnification  factors. 

The  fact  is  that  the  assumption  of  linear  response  is  virtually  never 
rigorously  justifiable  in  the  cases  of  greatest  practical  interest,  which 
are,  of  course,  the  cases  of  severe  motions.  Quite  the  contrary,  one 
usually  knows  very  well  that  the  responses  are  not  linear  in  cases  of 
large-amplitude  motions.  Nevertheless,  there  is  little  practical  recourse 
but  to  make  the  assumption  anyway.  The  real  justification  is  made  on 
the  practical  grounds  that  the  results  have  usually  been  found  to  agree 
sufficiently  well  with  experience  to  be  extremely  useful. 


Our  confidence  In  the  method,  however,  must  be  kept  somewhat  In 
proportion  to  the  extent  of  the  closely  related  successful  experience 
which  has  been  Accumulated.  The  experience  with  energy-spectral  anal- 
ysis of  Ground  Effect  Machine  motions  has  certainly  been  extremely 
limited,  at  best,  and  judgment  must  be  reserved. 

GEM's  with  significant  air  gaps  appear,  at  first  glance,  to  be 
quite  unsulted  to  this  type  of  analysis  because  of  the  obvious  major 
differences  between  the  character  of  small-amplitude  motions,  with  no 
water  contact,  and  motions  of  larger  amplitude  with  frequent  water  con- 
tact. On  the  other  hand,  in  severe  wave  conditions  the  water  contact 
may  be  nearly  continuous,  and  it  is  not  unlikely  that  useful  energy- 
spectral  analyses  of  such  motions  can  be  made,  if  appropriate  approxi- 
mations for  the  magnification  factors  can  be  established.  GEM's  of  the 
Captured  Air  Bubble  type,  or  other  types  which  operate  without  signifi- 
cant air  gaps,  would  seem  to  lend  themselves  very  readily  to  this  type 
of  analysis.  However,  this  also  remains  to  be  seen. 

Magnification  factors  are  most  frequently  determined  from  calcula- 
tions, from  towing-tank  tests  per#>rmed  in  regular  waves,  or  both.  The 
calculations  can  be  made  fairly  successfully  for  simple  displacement 
hull  forms,  but  the  towing  tank  is  relied  upon  for  less  conventional 
types  of  vehicles,  and  even  for  more  sophisticated  forms  of  displace- 
ment hulls. 

In  Section  F,  simplified  calculations  are  carried  out  for  the  heave 
magnification  factor  for  a full-peripheral  craft  with  no  water  contact. 

In  the  case  — > > 1 Equations  [F-3a]  and  [F-3bl  are  derived  and  then 

D 

a formula  is  derived  for  an  attenuation  factor  (Figure  F-3),  allowing 
for  effects  of  finite  wave  length.  Those  results  are  of  limited  practi- 
cal interest  becausei  (a)  the  accuracy  of  the  various  estimates  and 
assumptions  was  not  necessarily  very  good,  and  (b)  interest  in  "rough- 
water"  motions  without  water  contact  is  largely  academic,  anyway.  However, 
that  development  does,  at  least,  give  some  idea  of  what  is  involved  in 
calculating  the  magnification  factors.  It  is  a fairly  complex  problem, 
even  when  one  chosses  the  simplest  cases,  and  makes  very  liberal  use  of 
simplifying  approximations  It  is  frequently  next  to  impossible,  when 
one  comes  to  problems  of  real  practical  interest. 


r 


A typical  towlng-tank  determination  of  the  magnification  factors 
might  be  described  as  follows.  The  frequency  of  encounter  with  a simple 
wave  is 

Pf*  + 2n  jr 


Therefore,  at  given  model  velocity,  the  desired  range  of  frequency  of 
encounter  can  be  traced  out  by  testing  over  simple  waves  of  various 
wave  lengths,  thus  obtaining  curves  of  M^tu)  and  My(u))  at  various  con- 
stant velocities  V . These  must  then  be  adjusted  to  allow  for  the  scale 
o 

of  the  model,  In  accordance  with  Froude  scaling  laws.  The  magnification 
factors  are  nondimens ional , and  require  no  adjustment,  but  u>  and  VQ  must 
be  adjusted  In  a manner  which  can  be  seen  by  rewriting  the  above  expres- 
sion for  U)  In  nondlmens  ional  form: 


To  preserve  geometric  scale,  obviously 


model 


prototype 


Then,  since  the  acceleration  of  gravity  g is  Independent  of  scale, 
we  obtain 


prototype 


model 


(i) 

prototype 

model 


* 


<«») 


prototype 


^ model 


model 

( l ) 

v prototype 
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(The  nond linens lonal  velocity  parameter 


la  called  the  "Froude  number."  Since  g la  Invariant,  naval  architects 
frequently  use  "speed-length  ratio"  V^/*/T~  as  a more  convenient  velocity 
parameter . ) 

Having  adjusted  the  model  results  to  full  scale,  the  magnification 

factors  for  each  velocity  are  multiplied  by  the  appropriate  energy 

2 

spectrum  of  the  waves  [A(uu)  ] to  obtain  the  energy  spectra  of  the  motions, 
as  previously  discussed. 

The  necessary  hypothesis  that  the  responses  are  linear  Is  easily 
examined  during  the  model  tests,  by  determining  whether  the  motions  are 
sinusoidal,  and  whether  the  magnification  factors  are  Independent  of  wave 
height.  Of  course,  they  often  are  not,  with  any  very  high  degree  of 
exactness.  The  questions  to  be  answered  then  are:  (a)  whether  the  non- 
linearities  are  sufficiently  sligtit  that  one  can  proceed  with  energy- 
spectral  analysis,  and  (b)  If  so,  how  should  one  choose  the  best  "equiv- 
alent linear"  approximations  to  the  magnification  factors,  from  the 
experimental  measurements  obtained.  These  questions  are  best  left  to  be 
answered,  on  a case-by-case  basis,  by  veteran  experimenters. 

The  preceding  discussion  has  been  confined  to  prediction  of  the 
pitch  and  heave  energy  spectra  for  vehicles  In  one-dimensional  head  or 
following  seas.  The  extension  to  prediction  of  other  properties  of  the 
motion  is  obvious,  so  long  as  it  is  possible  to  determine  linear-response 
magnification  factors  for  the  desired  properties.  ♦ 

The  extensions  to  motions  at  headings  other  than  normal  to  wave 
crests,  and  to  motions  over  two-dimensional  wave  systems  (i.e.,  systems 
composed  of  simple-wave  components  propagating  in  a variety  of  directions) 
are  beyond  the  scope  of  our  Introductory  discussion. 

VA-3  TESTS 

Considering  the  Importance  which  is  attached  to  seakeeping  by  those 
engaged  in  GEM  research  and  development,  there  Is  a remarkable  paucity 
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of  experimental  information  on  this  subject  available  in  open  literature. 
To  gome  extent,  this  reflects  the  fact  that  a thorough  evaluation  of  the 
seakeeping  characteristics  of  even  a single  design  is  a very  difficult, 
expensive,  and  time-consuming  task.  More  importantly,  it  undoubtedly  re- 
flects the  fact  that  the  more  important  investigations  have  been  slow  in 
finding  their  way  into  the  literature,  for  commercial  reasons. 

The  most  important  data  available  in  open  literature  nj:  present 
appear  to  be  the  results  of  tests  of  the  Vickers  VA-3  craft  performed 
by  Republic  Aviation  Corporation,  Reference  2.  These  tests  were  quite 
limited  in  scope,  but  sufficient  results  were  obtained  to  give  some  idea 
of  the  riding  qualities  of  a peripheral-jet  type  machine  with  flexible 
trunks  and  small  daylight  gap. 

The  results,  in  addition  to  the  quantitative  information  they  pro- 
vide, are  of  great  value  in  confirming  numerous  favorable,  but  largely 
unsupported  statements  in  the  press  regarding  Hovercraft  seakeeping 
qualities. 

A general-arrangement  drawing  of  the  VA-3  is  reproduced  from  Refer- 
ence 2 (see  Figure  G-l).  The  VA-3  had  a gross  weight  of  about  30,000 
pounds,  cushion  area  of  about  1000  square  feet,  and  mean-cushion-length/ 
beam  ratio  l/b  of  about  2.0.  (These  dimensions  are  roughly  estimated 
from  the  drawing.) 

A representative  set  of  statistical  data  on  vertical  accelerations 
measured  at  the  bow  are  reproduced  in  Figure  G-2.  Data  allowing  compar- 
ison of  the  acceleration  levels  at  the  bow  with  those  at  the  center  of 
gravity  are  reproduced  in  Figure  G-3. 

These  data  pertain  to  the  maxima  of  the  acceleration  time-histories, 
the  ordinates  F of  the  curves  giving  the  probability  that  an  arbitrarily 
chosen  maximum  would  exceed  the  corresponding  abscissas.  Thus  the 
abscissa  corresponding  to  F ■ 0.5  is  the  median  of  all  the  maxima;  the 
abscissa  corresponding  to  P » 7(7)  ° *s  t^le  median  of  the  one- 

third.  highest  maxima,  etc.  Reference  2 states  that  the  statistical 


quality  of  the  data  la  not  very  high,  due  to  comparatively  rapid  changes 
in  the  energy  spectrum  of  the  waves  and  imprecise  control  of  other  test 


w 
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variables.  This  Is  reflected  In  the  Irregular  shapes  of  the  probability 
curves  and  some  seeming  Inconsistencies  between  results  for  different 
conditions.  Nevertheless,  the  results  provide  a very  useful  Indication 
of  the  severity  of  motions  likely  to  be  encountered. 

The  results  can  be  rendered  somewhat  more  useful  by  representing 
them  In  an  appropriate  nondimens ional  form.  This  Is  done  In  Figure  G-4, 
where  the  entire  envelope  of  the  VA-3  tests  is  shown  on  a map  of  speed- 
length  ratio  y/J  Z versus  ratio  of  mean  wave  height  to  cushion  length 

H / F , and  then  the  "maximum"  bow  accelerations  (median  of  the  one-tenth 

av 

highest  maxima)  for  various  conditions  represented  In  Figure  G-2  are 
located  on  this  map.  (As  previously  noted,  the  speed-length  ratio  V^/ J 1 

can  be  considered  nondlmens ional  If  the  acceleration  of  gravity  g Is 
considered  constant.)  Obviously,  the  results  are  Insufficient  to  permit 
any  detailed  trends  to  be  mapped.  The  results  are  significant,  however, 
as  a qualitative  Indication  of  the  remarkable  ability  of  GEM's  to  traverse 
rough  water  at  speed  without  severe  motions. 

These  results  are  particularly  Impressive  when  the  nondlmens ional 
map  of  Figure  G-4  is  re-interpreted  in  dimensional  form  for  vehicles 
larger  than  VA-3.  (This  Is  strictly  valid  provided  the  larger  vehicles 
are  "scale  models"  of  VA-3  and  the  scaled-up  sea  has  proper  spectrum,  In 
accordance  with  the  Froude  scaling  rules.)  An  extreme  example  is  given  on 
Figure  G-2  In  the  form  of  auxiliary  scales  along  the  top  and  right-hand 
side  of  the  graph  giving  wave  height  and  speed  for  a "model"  8.9  times  as 
large  as  VA-3,  weighing  (8.9)3  times  as  much.  Keep  In  mind  that  the 
accelerations  mapped  pertain  to  the  bow,  the  accelerations  at  the  c.g. 
being  only  about  40  percent  as  great. 

Of  course,  when  9400-ton  GEM's  are  built,  they  are  quite  unlikely  to 
be,  strictly  speaking,  "scale  models"  of  VA-3.  Nevertheless,  one  Is  strongly 
encouraged  to  conclude  that  seakeeping  problems  would  not  constitute  a 
serious  barrier  to  the  development  and  application  of  large  GEM's. 

Aerodynamics  Laboratory 
David  Taylor  Model  Basin 
Washington,  D.  C.  20007 
July  1966 
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Figure  G-l  - General  Arrangement  of  VA-3 


Vertical  Accelerations  of  VA-3  at 
(a)  Sea  State  2 


(F  - Fraction  of  acceleration  peaks  exceeding  given  value.) 

Figure  G-2  (Concluded) 

(b)  Sea  State  3 


(F  - Fraction  of  acceleration  peaks  exceeding  given  value.) 

Figure  G-3  - Comparison  of  Acceleration  at  Bow  and  Center  of  Gravity 


Speed-Length  Ratio 
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Wave  Height  Ratio, 


Figure  G-4  - Approximate  Map  of  VA-3  Teat  Results,  h/i  = 0.005 

(Reference  2) 


26 


in  knots  for  9,400-ton  "Scale  Model"  of  VA-3 
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